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Abstract 

, We present the theory, the experimental evidence and fundamental phys- 

ical consequences concerning the existence of families of undistorted progres- 
sive waves (UPWs) of arbitrary speeds < v < oo, which are solutions of 
\£) . the homogeneous wave equation, Maxwell equations, Dirac, Weyl and Klein- 

Gordon equations. 



PACS numbers: 41.10.Hv; 03.30.+p; 03.40Kf 
1. Introduction 



In this paper we present the theory, the experimental evidence, and the fun- 
damental physical consequences concerning the existence of families of undistorted 
progressive waves (UPWs)*-*-' moving with arbitrary speeds < -**- 1 < v < oo. We show 
that the main equations of theoretical physics, namely: the scalar homogeneous 
wave equation (HWE); the Klein-Gordon equation (KGE); the Maxwell equations, 
the Dirac and Weyl equations have UPWs solutions in a homogeneous medium, in- 
cluding the vacuum. By UPW, following Courant and Hubert' 1 ' we mean that the 
UPW waves are distortion free, i.e. they are translationally invariant and thus do 
not spread, or they reconstruct their original form after a certain period of time. 
Explicit examples of how to construct the UPWs solutions for the HWE are found in 
Appendix A. The UPWs solutions to any field equations have infinite energy. How- 
ever, using the finite aperture approximation (FAA) for diffraction (Appendix A), 

(*)tjPW is used for the singular, i.e., for undistorted progressive wave. 
(**)'We use units where c = 1, c being the so called velocity of light in vacuum. 
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we can project quasi undistorted progressive waves (QUPWs) for any field equation 
which have finite energy and can then in principle be launched in physical space. 

In section 2 we show results of a recent experiment proposed and realized by 
us where the measurement of the speeds of a FAA to a subluminal^ Bessel pulse 
[eq.(2.1)] and of the FAA to a superluminal X -wave [eq.(2.5)] are done. The results 
are in excellent agreement with the theory. 

In section 3 we discuss some examples of UPWs solutions of Maxwell equa- 
tions; (i) subluminal solutions which are interesting concerning some recent attempts 
appearing in the literature' 2 ' 3 ' 4 ^ of construction of purely electromagnetic particles 
(PEP) and (ii) a superluminal UPW solution of Maxwell equations called the su- 
perluminal electromagnetic X-wave^ (SEXW). We briefly discuss how to launch a 
FAA to SEXW. In view of the experimental results presented in section 2 we are 
confident that such electromagnetic waves will be produced in the next few years. In 
section 4 we discuss the important question concerning the speed of propagation of 
the energy carried by superluminal UPWs solutions of Maxwell equations, clearing 
some misconceptions found in the literature. In section 5 we show that the experi- 
mental production of a superluminal electromagnetic wave implies in a breakdown 
of the Principle of Relativity. In section 6 we present our conclusions. 

Appendix B presents a unified theory of how to construct UPWs of arbitrary 
speeds < v < oo which are solutions of Maxwell, Dirac and Weyl equations. Our 
unified theory is based on the Clifford bundle formalism' 6 ' 7 ' 8 ' 9 ' 10 ^ where all fields 
quoted above are represented by objects of the same mathematical nature. We take 
the care of translating all results in the standard mathematical formalisms used by 
physicists in order for our work to be usefull for a larger audience. 

Before starting the technical discussions it is worth to briefly recall the history of 
the UPWs of arbitrary speeds < v < oo, which are solutions of the main equations 
of theoretical physics. 

To the best of our knowledge H. Bateman' 11 ' in 1913 was the first person to 
present a subluminal UPW solution of the HWE. This solution corresponds to what 
we called the subluminal spherical Bessel beam in Appendix A [see eq.(A.31)]. Ap- 
parently this solution has been rediscovered and used in diverse contexts many 
times in the literature. It appears, e.g., in the papers of Mackinnon' 12 ! of 1978 and of 
Gueret and Vigier' 13 ' and more recently in the papers of Barut and collaborators' 14 ' 15 ^ . 
In particular in' 14 ! Barut also shows that the HWE has superluminal solutions. In 

Win this experiment the waves are sound waves in water and, of course, the meaning of the 
words subluminal, luminal and superluminal in this case is that the waves travel with speed less, 
equal or greater than c s , the so called velocity of sound in water. 
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1987 Durnin and collaborators rediscovered a subluminal UPW solution of the HWE 
in cylindrical coordinates' 16,17,18 '. These are the Bessel beams of section A4 [see 
eq.(A.41)]. We said rediscovered because these solutions are known at least since 
1941, as they are explicitly written down in Stratton's book' 19 '. The important point 
here is that Durnin' 16 ' and collaborators constructed an optical subluminal Bessel 
beam. At that time they didn't have the idea of measuring the speed of the beams, 
since they were interested in the fact that the FAA to these beams were quasi UPWs 
and could be very usefull for optical devices. Indeed they used the term "diffraction- 
free beams" which has been adopted by some other authors later. Other authors 
still use for UPWs the term non- dispersive beams. We quote also that Hsu and 
collaborators' 20 ' realized a FAA to the Jo Bessel beam [eq.(A.41)] with a narrow band 
PZT ultrasonic transducer of non-uniform poling. Lu and Greenleaf' 21 ' produced the 
first Jo nondiffr acting annular array transducers with PZT ceramic/polymer com- 
posite and applied it to medical acoustic imaging and tissue characterization' 22,23 ' . 
Also Campbell et al' 24 ' used an annular array to realize a FAA to a J Bessel beam 
and compared the Jo beam to the so called axicon beam^ 25 \ For more on this topic 
seel 26 '. 

Luminal solutions of a new kind for the HWE and Maxwell equations, also 
known as focus wave mode [FWM] (see Appendix A), have been discovered by 
Brittingham ' 27 ' (1983) and his work inspired many interesting and important studies 
as, e.#.,[ 29 - 4 °]. 

To our knowledge the first person to write about the possibility of a superluminal 
UPW solution of HWE and, more important, of Maxwell equations was Band' 41,42 '. 
He constructed a superluminal electromagnetic UPW from the modified Bessel beam 
[eq.(A.42)j which was used to generate in an appropriate way an electromagnetic 
potential in the Lorentz gauge. He suggested that his solution could be used to 
eventually launch a superluminal wave in the exterior of a conductor with cylin- 
drical symmetry with appropriate charge density. We discuss more some of Band's 
statements in section 4. 

In 1992 Lu and Greenleaf' 43 ' presented the first superluminal UPW solution 
of the HWE for acoustic waves which could be launched by a physical device' 44 '. 
They discovered the so called X- waves, a name due to their shape (see Fig. 3). 
In the same year Donnelly and Ziolkowski' 45 ' presented a thoughtfull method for 
generating UPWs solutions of homogeneous partial equations. In particular they 
studied also UPW solutions for the wave equation in a lossy infinite medium and to 
the KGE. They clearly stated also how to use these solutions to obtain through the 
Hertz potential method (see appendix B, section B3) UPWs solutions of Maxwell 
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equations. 

In 1993 Donnely and Ziolkowski' 46 ' reinterpreted their study off 45 ' and obtained 
subluminal, luminal and superluminal UPWs solutions of the HWE and of the KGE. 
In Appendix A we make use of the methods of this important paper in order to obtain 
some UPWs solutions. Also in 1992 Barut and Chandola' 47 ! found superluminal 
UPWs solutions of the HWE. In 1995 Rodrigues and VW 48 ! discovered in quite 
an independent way^ subluminal and superluminal UPWs solutions of Maxwell 
equations and the Weyl equation. At that time Lu and Greenleaf® proposed also 
to launch a superluminal electromagnetic X-wave/**- 1 

In September 1995 Professor Ziolkowski took knowledge of' 48 ^ and informed one 
of the authors [WAR] of his publications and also of Lu's contributions. Soon a 
collaboration with Lu started which produced this paper. To end this introduction 
we must call to the reader's attention that in the last few years several important ex- 
periments concerning the superluminal tunneling of electromagnetic waves appeared 
in the literature^ 51 ' 52 ]. Particularly interesting is Nimtz's paper' 53 ! announcing that 
he transmitted Mozart's Symphony # 40 at 4.7c through a retangular waveguide. 
The solutions of Maxwell equations in a waveguide lead to solutions of Maxwell 
equations that propagate with subluminal or superluminal speeds. These solutions 
can be obtained with the methods discussed in this paper and will be discussed in 
another publication. 

2. Experimental Determination of the Speeds of Acoustic 
Finite Aperture Bessel Pulses and X- Waves. 

In appendix A we show the existence of several UPWs solutions to the HWE, 
in particular the subluminal UPWs Bessel beams [eq.(A.36)j and the superluminal 
UPWs X-waves [eq.(A.52)j. Theoretically the UPWs X-waves, both the broad-band 
and band limited [see eq.(2.4)] travel with speed v = c s / cost] > 1. Since only FAA 
to these X-waves can be launched with appropriate devices, the question arises if 
these FAA X-waves travel also with speed greater than c s , what can be answered 
only by experiment. Here we present the results of measurements of the speeds of a 

(*)Rodrigues and Vaz are interested in obtaining solutions of Maxwell equations characterized 
by non-null field invariants, since solutions of this kind are! 49,50 ' necessary in proving a surprising 
relationship between Maxwell and Dirac equations. 

(**) A 

version of [5] was submitted to IEEE Trans. Antennas Propag. in 1991. See reference 

40of[ 43 l 
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FAA to a broad band Bessel beam, called a Bessel pulse (see below) and of a FAA 
to a band limited X-wave, both moving in water. We write the formulas for these 
beams inserting into the HWE the parameter c s known as the speed of sound in 
water. In this way the dispersion relation [eq.(A.37)j must read 

^-k 2 = a\ (2.1) 

Then we write for the Bessel beams 

$<(*,£) = Uap)e^- t+ne \ n = 0,l,2,... (2.2) 

Bessel pulses are obtained from eq.(2.2) by weighting it with a transmitting transfer 
function, T{uj) and then linearly superposing the result over angular frequency uj, 
i.e., we have 

^j BBn (t,x) = 27ce me .Uap)F~ 1 iT(uy k % (2.3) 

where T~ x is the inverse Fourier transform. The FAA to &j BB will be denoted by 
FAA$< BBn (ord>< AJ J. 

We recall that the X-waves are given by eq.(A.52), i.e., 

POO — 

<S>> n (t,x)=e ine B(k)J n (kpsmr])e- k[aa - l{zcosv - Cst)] dk, (2.4) 

J 

where k — kj cosr/, k = uo/c s . By choosing B(k) = a we have the infinite aperture 
broad bandwidth X-wave [eq.(A.53)j given by 



M(r + y/M) n ' 

(2.5) 

M = (psmri) 2 + r 2 , r = [a — i(zcosi] — c s t)]. 

A FAA to §xBB n wil1 be denoted by FAA§^ BBn . When B(k) in eq.(2.4) is different 
from a constant, e.g., if B(k) is the Blackman window function we denote the X- 
wave by $x_BL n > where BL means band limited. A FAA to &xBL n wn l be denoted 
FAA$x_BL n - Also when T(uo) in eq.(2.3) is the Blackman window function we denote 
the respective wave by §jbl„- 

As discussed in Appendix A and detailed in' 26 ' 44 ] to produce a FAA to a given 
beam the aperture of the transducer used must be finite. In this case the beams 
produced, in our case FAA& JBLo and FAA$xbb , have a finite depth of field l 26 ' 
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(DF)W and can be approximately produced by truncating the infinite aperture 
beams Qjbl and &xbb (or &xbl ) at the transducer surface (z — 0). Broad band 
pulses for z > can be obtained by first calculating the fields at all frequencies with 
eq.(A.28), i.e., 



~ | ra ~ gikR 

<& FA (u,x) = —J J p'dp'd6'$(uj,x')—z 

(2.6) 

where the aperture weighting function $(u;, x') is obtained from the temporal Fourier 
transform of eqs.(2.3) and (2.4). If the aperture is circular of radius a [as in eq.(2.6)], 
the depth of field of the FAA$ jbl pulse, denoted BZ max and the depth of field of 
the FAA$xbb or FAA $ X bl denoted by XZ max are given by^ 261 




BZ max = a W - 1 ; XZ max = a cot rj. (2.7) 



For the FAA<3>j Bio pulse we choose T(uo) as the Blackman window function^ 54 ! 
that is peaked at the central frequency f = 2.5MHz with a relative bandwidth of 
about 81% (—6dB bandwidth divided by the central frequency). We have 



B(k) = < 



a 



0.42 -0.5^ + 0.08 cos 2:VA 



k k 
otherwise. 



, < k < 2k ; (2 g) 



The "scaling factor" in the experiment is a = 1202.45m- 1 and the weighting function 

&jbb (u, x) in eq.(2.6) is approximated with stepwise functions. Practically this is 
done with the 10-element annular array transfer built by Lu and Greenleaf' 26 ' 44 l 
The diameter of the array is 50mm. Fig. 1***) shows the block diagram for the 
production of FAA &xbl an d FAA$j BLq . The measurement of the speed of the 
FAA Bessel pulse has been done by comparing the speed with which the peak of the 
FAA Bessel pulse travels with the speed of the peak of a pulse produced by a small 



(*)DF is the distance where the field maximum drops to half the value at the surface of the 
transducer. 

("^Reprinted with permission from fig. 2 of! 44 '. 
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circular element of the array (about 4mm or 6.67A in diameter, where A is 0.6mm 
in water). This pulse travels with speed c s = 1.5mm/ fis. The distance between 
the peaks and the surface of the transducer are 104.33(9)mm and 103.70(5)mm for 
the single-element wave and the Bessel pulse, respectively, at the same instant t of 
measurement. The results can be seen in the pictures taken from of the experiment 
in Fig. 2. As predicted by the theory developed in Appendix A the speed of the 
Bessel pulse is 0.611(3)% slower than the speed c s of the usual sound wave produced 
by the single element. 

The measurement of the speed of the central peak of the FAA Q X bl wave 
obtained from eq.(2.4) with a Blackman window function [eq.(2.8)] has been done in 
the same way as for the Bessel pulse. The FAA$xrl wave has been produced by 
the 10-element array transducer of 50mm of diameter with the techniques developed 
by Lu and Greenleaf' 26 ' 44 ! . The distances traveled at the same instant t by the single 
element wave and the X-wave are respectively 173.48(9)mm and 173.77(3)mm. Fig. 
3 shows the pictures taken from the experiment. In this experiment the axicon angle 
is i] = 4°. The theoretical speed of the infinite aperture X-wave is predicted to be 
0.2242% greater then c s . We found that the FAA$xbb wave traveled with speed 
0.267(6)% greater then c s ! 

These results, which we believe are the first experimental determination of the 
speeds of subluminal and superluminal quasi-UPWs FAA$ j BL and FAA$ j BBo solu- 
tions of the HWE, together with the fact that, as already quoted, Dunlin^ 16 ' produced 
subluminal optical Bessel beams, give us confidence that electromagnetic subluminal 
and superluminal waves may be physically launched with appropriate devices. In 
the next section we study in particular the superluminal electromagnetic X-wave 



It is important to observe here the following crucial points: (i) The FAA &xbb„ is 
produced by the source (transducer) in a short period of time At. However, different 
parts of the transducer are activated at different times, from to At, calculated from 
eqs.(A.9) and (A. 28). As a result the wave is born as an integral object for time 
At and propagates with the same speed as the peak. This is exactly what has been 
seen in the experiments and is corroborated by the computer simulations we did for 
the superluminal electromagnetic waves (see section 3). (ii) One can find in almost 
all textbooks that the velocity of transport of energy for waves obeying the scalar 
wave equation 



(SEXW). 




(2.9) 
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(2.10) 



(2.11) 



from which it follows that 




(2.12) 



Our acoustic experiment shows that for the X-waves the speed of transport of 
energy is c s / cost], since it is the energy of the wave that activates the detector 
(hydro-phone). This shows explicitly that the definition of v £ is meaningless. This 
fundamental experimental result must be kept in mind when we discuss the meaning 
of the velocity of transport of electromagnetic waves in section 4. 



Figure 1: Block diagram of acoustic production of Bessel pulse and X- Waves. 
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Figure 2: Propagation speed of the peak of Bessel pulse and its comparison with that 
of a pulse produced by a small circular element (about 4 mm or 6.67 A in diameter, 
where A is 0.6 mm in water). The Bessel pulse was produced by a 50 mm diameter 
transducer. The distances between the peaks and the surface of the transducer 
are 104.339 mm and 103.705 mm for the single-element wave and the Bessel pulse, 
respectively. The time used by these pulses is the same. Therefore, the speed of the 
peak of the Bessel pulse is 0.611(3)% slower than that of the single-element wave. 
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Figure 3: Propagation speed of peak of X-wave and its comparison with that of 
a pulse produced by small circular element (about 4 mm or 6.67 A, where A is 
0.6 mm in water). The X-wave was produced by a 50 mm diameter transducer. 
The distance between the peaks and the surface of the transducer are 173.489 mm 
and 173.773 mm for the single-element wave and the X-wave, respectively. The time 
used by these pulses is the same. Therefore, the speed of the peak of the X-wave 
is 0.2441(8)% faster than that of the single-element wave. The theoretical ratio for 

( C / COS 71 — C ) 

X-waves and the speed of sound is — - = 0.2442% for rj = 4°. 
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3. Subluminal and Superluminal UPWs Solutions of Maxwell 
Equations(ME) 

In this section we make full use of the Clifford bundle formalism (CBF) resumed 
in Appendix B, but we give translation of all the main results in the standard vec- 
tor formalism used by physicists. We start by reanalyzing in section 3.1 the plane 
wave solutions (PWS) of ME with the CBF. We clarify some misconceptions and ex- 
plain the fundamental role of the duality operator 75 and the meaning of i = y/—T in 
standard formulations of electromagnetic theory. Next in section 3.2 we discuss sub- 
luminal UPWs solutions of ME and an unexpected relation between these solutions 
and the possible existence of purely electromagnetic particles (PEPs) envisaged by 
Einstein t 55 ', Poincare^ 56 ', Ehrenfestt 57 ] and recently discussed by Waite, Barut and 
Zeni' 2,3 !. In section 3.3 we discuss in detail the theory of superluminal electromag- 
netic X-waves (SEXWs) and how to produce these waves by appropriate physical 
devices. 

3.1 Plane Wave Solutions of Maxwell Equations 

We recall that Maxwell equations in vacuum can be written as [eq.(B.6)j 

<9F = 0, (3.1) 

where Fsec/\ 2 (M) C sec C£(M). The well known PWS of eq.( |3.1| ) are obtained 
as follows. We write in a given Lorentzian chart of the maximal atlas of M 
(section B2) a PWS moving in the z-direction 

F = fe< hkx , (3.2) 

k = ^7 M , k 1 = k 2 = 0, x = x^, (3.3) 

where k, x G secA^-^O C sec C£(M) and where / is a constant 2-form. From 
eqs. (|3TID and ( |3.2| ) we obtain 

kF = (3.4) 

Multiplying eq.(3.4) by k we get 

k 2 F = (3.5) 
and since k G secA 1 (M) C sec C£(M) then 

k 2 = <-> k° = ±\k\ = ±k 3 , (3.6) 
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i.e., the propagation vector is light-like. Also 

F 2 = F. F + F AF = (3.7) 

as can be easily seen by multiplying both members of eq.(3.4) by F and taking into 
account that k 7^ 0. Eq(3.7) says that the field invariants are null. 

It is interesting to understand the fundamental role of the volume element 75 
(duality operator) in electromagnetic theory. In particular since e l5kx = cos kx + 

— * — * 

75 sin kx, 75 = i, writing F = E + '\B (see eq.(B.17)), / = e± + ie 2 , we see that 

E + \B = e\ cos kx — e 2 sin kx + \{e\ sin kx + e 2 cos kx) . (3.8) 

— * — * 

From this equation, using OF = 0, it follows that e*i.e 2 = 0, k.t\ = k.e 2 = and 
then 

E.B = 0. (3.9) 

This equation is important because it shows that we must take care with the i = 
\J — 1 that appears in usual formulations of Maxwell Theory using complex electric 
and magnetic fields. The % = yf^l in many cases unfolds a secret that can only be 
known through eq.(3.8). It also follows that k.E = k.B = 0, i.e., PWS of ME are 
transverse waves. We can rewrite eq.(3.4) as 

^7070^70 = (3.10) 

and since k^o = ko + k, 70-^70 = — E + \B we have 

kf = k f. (3.11) 

Now, we recall that in C£ + (M) (where, as we say in Appendix B, the typical 
fiber is isomorphic to the Pauli algebra C£ 3i0 ) we can introduce the operator of space 

— * 

conjugation denoted by * such that writing / = e + ib we have 

/* = -e + i6 ; k* = k ■ k* = -k. (3.12) 

We can now interpret the two solutions of k = 0, i.e. k = \k\ and k = — |fc| as 

— * — * 

corresponding to the solutions k f = kf and k f* = —kf*; f and /* correspond 
in quantum theory to "photons" of positive or negative helicities. We can interpret 
ko — \k\ as a particle and ko = — \k\ as an antiparticle. 

Summarizing we have the following important facts concerning PWS of ME: (i) 
the propagation vector is light-like, k 2 = 0; (ii) the field invariants are null, F 2 = 0; 
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(iii) the PWS are transverse waves, i.e., k.E = k.B = 0. 

3.2 Subluminal Solutions of Maxwell Equations and Purely Electromag- 
netic Particles. 

We take $ £ sec(A°(M) © A 4 (^0) C sec Cl{M) and consider the following Hertz 
potential vr £ secA 2 (M) C secC£(M) [eq. flggg )] 

7r = $ 7 1 7 2 . (3.13) 

We now write 

$(t,x) = 0(f)e 7sOi . (3.14) 
Since 7r satisfies the wave equation, we have 

V 2 (f)(x) +n 2 (j)(x) = 0. (3.15) 

Solutions of eq.(3.15) (the Helmholtz equation) are well known. Here we consider 
the simplest solution in spherical coordinates, 

0(f) = C^^- , r = ^f x 2 + y 2 + z 2 , (3.16) 

where C is an arbitrary real constant. From the results of Appendix B we obtain 
the following stationary electromagnetic field, which is at rest in the reference frame 
Z where are naturally adapted coordinates (section B2): 

C 

F = — [sin Qt(aQr sin 8 simp — j3 sin 8 cos 8 cos ^'yoli 
— sin Qt(aflr sin 8 cos </? + (3 sin 8 cos 8 sin 95)7072 
+ sin Qt(P sin 2 8 - 2a)7 7 3 + cos flt((3 sin 2 8 - 2a)^ 2 (3. 17) 
+ cos Qt(j3 sin 8 cos 8 sin ip + aflr sin 8 cos 92)7173 
+ cos Qt(—(3 sin 8 cos 8 cos <p + aflr sin 8 sin 92)7273] 

with a = fir cos fir — sin fir and f3 = 3a + Q 2 r 2 sin fir. Observe that Fq is regular at 
the origin and vanishes at infinity. Let us rewrite the solution using the Pauli-algebra 
in C£ + (M). Writing (i = 75) 

F = E Q + iB (3.18) 

we get 
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E = W sin Vlt, B = W cos fit, (3. 19) 

with 

W = -C ( — ^ - — - — P^ + y 2 } _ (3 20) 

We verify that divW> = 0, div£ = di Vj B = 0, rotE +dB /dt = 0, rotB -dE /dt = 
0, and 

rotVV = fiW>. (3.21) 

Now, from eq.(B.88) we know that T = -F'-foF is the 1-form representing the 

energy density and the Poynting vector. It follows that EqxBq = 0, i. e., the solution 
has zero angular momentum. The energy density u = S 00 is given by 



u = ^ [sin 2 fl(fi 2 rV + (3 2 cos 2 9) + ((3 sin 2 9 - 2a) 2 ] (3.22) 



Then J J j R3 udv — oo. As explained in section A. 6 a finite energy solution can be 
constructed by considering "wave packets" with a distribution of intrinsic frequencies 
F(fi) satisfying appropriate conditions. Many possibilities exist, but they will not 
be discussed here. Instead, we prefer to direct our attention to eq.(3.21). As it 
is well known, this is a very important equation (called the force free equation' 2 ') 
that appears e.g. in hydrodynamics and in several different situations in plasma 
physics' 58 '. The following considerations are more important. 

Einstein' 55 ! among others (see' 3 ' for a review) studied the possibility of construct- 
ing purely electromagnetic particles (PEPs). He started from Maxwell equations for 
a PEP configuration described by an electromagnetic field F p and a current density 
J p , where 

dF p = J p (3.23) 
and rightly concluded that the condition for existence of PEPs is 

Jp.Fp = 0. (3.24) 

This condition implies in vector notation 

PpE p = 0, jp.E p = 0, j p xBp = 0. (3.25) 
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From eq.(3.24) Einstein concluded that the only possible solution of eq.(3.22) with 
the subsidiary condition given by eq.(3.23) is J p = 0. However, this conclusion is 
correct, as pointed in' 2 ' 3 !, on iy jf j2 > q ; i.e., if J p is a time-like current density. 
However, if we suppose that J p can be spacelike, i.e., J 2 p < 0, there exists a reference 
frame where p p = and a possible solution of eq.(3.24) is 

Pp = 0, E p .B p = 0, j p = KCB p , (3.26) 
where K = ±1 is called the chirality of the solution and C is a real constant. In' 2,3 ] 

— * 

static solutions of eqs.(3.22) and (3.23) are exhibited where E p = 0. In this case we 

— * 

can verify that B p satisfies 

V x B p = KCB P . (3.27) 
Now, if we choose F e sec/\ 2 (M) C sec COM) such that 

E Q = B p cos Qt, Bq = Bp sin Qt 

and fl = KC > 0, we immediately realize that 

dF = 0. (3.29) 

This is an amazing result, since it means that the free Maxwell equations may 
have stationary solutions that may be used to model PEPs. In such solutions the 
structure of the field F is such that we can write 

F = F' + F = iWcosnt-Wsinnt, , 

OFp = -dF = J'p, (3 - 30) 

i. e., dF = is equivalent to a field plus a current. This fact opens several interesting 
possibilities for modeling PEPs (see also I 4 ') and we discuss more this issue in another 
publication. 

We observe that moving subluminal solutions of ME can be easily obtained 
choosing as Hertz potential, e.g., 

sin ^~2^<^ 

n<{t,x) = C — - — ^exp[7 5 (w<t - k^)]^^, (3.31) 
col-k 2 . = fi 2 ; 

U = [x 2 + y 2 + l 2 < (z-v < t) 2 }, (3.32) 
7< = I v < = du < /dk < . 
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We are not going to write explicitly the expression for F K corresponding to 7r < 
because it is very long and will not be used in what follows. 

We end this section with the following observations: (i) In general for sublu- 
minal solutions of ME (SSME) the propagation vector satisfies an equation like 
eq.(3.30). (ii) As can be easily verified, for a SSME the field invariants are non- 
null, (iii) A SSME is not a transverse wave. This can be seen explicitly from 
eq.(3.21). Conditions (i), (ii) and (iii) are in contrast with the case of the PWS of 
ME. W 49 ' 50 ! Rodri gues and Vaz showed that for free electromagnetic fields (OF = 0) 
such that F 2 ^ 0, there exists a Dirac-Hestenes equation (see section A. 8) for 
G sec(A°(M) + A 2 (M) + A 4 (M)) C sec C£(M) where F = ^7i72V>- This was the 
reason why Rodrigues and Vaz discovered subluminal and superluminal solutions of 
Maxwell equations (and also of Weyl equation) ' 48 J which solve the Dirac-Hestenes 
equation |eq.( [B.40|) l. 

3.3 The Superluminal Electromagnetic A- Wave (SEXW) 

To simplify the matter in what follows we now suppose that the functions 
[eq.(A.52)j and $xbb„ [eq.(A.53)j which are superluminal solutions of the scalar 
wave equation are 0- forms sections of the complexified Clifford bundle C£c{M) = 
(S ® Ct(M) (see section B4). We rewrite eqs.(A.52) and (A.53) as w 

/•oo _ 

$ x (t,x)=e ine B(k)J n (kpsmr])e- k[ao - i{zcos,1 - t)] dk (3.33) 
Jo 

and choosing B(k) = a , we have 

- ( . a (p sin v ) n e me 

ixBBjt ' x) = vmrrvw ( ' 

M = (psmrj) 2 + r 2 ; r — [ao — i(zcosr) — t)]. (3.35) 

As in section 2, when a finite broadband X-wave is obtained from eq.(3.31) with 
B{k) given by the Blackman spectral function [eq.(2.8)] we denote the resulting X- 
wave by &xbl„ {BL means band limited wave). The finite aperture approximation 
(FA A) obtained with eq.(A.28) to will be denoted FAA$xbl„ and the FAA 

to &xBB n will be denoted by FAA$ xbb„- We use the same nomenclature for the 
electromagnetic fields derived from these functions. Further, we suppose now that 



Win what follows n = 0, 1, 2, 
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the Hertz potential ir, the vector potential A and the corresponding electromagnetic 
field F are appropriate sections oiC£c(M). We take 

7T = $ 7l72 g sec (F ® A 2 (M) C sec C£ C {M), (3.36) 

where $ can be $ Xn , $jtBB„ , $xbl„ , FAA $xbb„ or FAA$ X bl„. Let us start by 
giving the explicit form of the F X bb„, i-e., the SEXWs. In this case eq.(B.81) gives 
n = n m and 

*m = ®XBB n Z (3.37) 

where z is the versor of the z-axis. Also, let p, be respectively the versors of the 
p and 9 directions where (p, 9, z) are the usual cylindrical coordinates. Writing 

— * — * 

FxBB n = E X BB n + lbB X BB n (3.38) 

we obtain from equations (A. 53) and (B.25): 

p d 2 d 2 
ExBB n = ---^*XBB n + ^$X BB „ 5 (3.39) 

d 2 Id 2 ( d 2 d 2 \ 

BxBBn = PQ^fxBB n + B- — * XBBn + Z \^*XBB n - ^XBB n j J (3.40) 

Explicitly we get for the components in cylindrical coordinates: 

(ExbbJ p = --u-^LQxbbJ (3.41a) 

(ExbbJo = -i^=—®xBB n ; (3.41b) 
p VMM 2 

(BxbbJp = cosr)(E X BB n )e] (3.41c) 

(BxBB n )e = ~cosr](E XBBn ) p ; (3.41d) 

(BxbbJ z = -sin 2 r 1 ^L§xBB n - (3.41e) 



The functions Mj, (i = 2, . . . , 7) in (3.41) are: 

M 2 = r + >/M; (3.42a) 

M 3 = n + -^r; (3.42b) 

M 4 = 2n + -jLr; (3.42c) 
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M 



M 5 = t + nVM; (3.42d) 
M 6 = (p 2 sin 2 - nM 3 )M 2 + np 2 ^- sin 2 77; (3.42e) 

M 7 = (n 2 - l)^ + 3n4Tr + 3^^r 2 . (3.42f) 

We immediately see from eqs.(3.41) that the F X bb„ are indeed superluminal 
UPWs solutions of ME, propagating with speed 1/ cos r] in the z-direction. That 
FxBB n are UPWs is trivial and that they propagate with speed c± — 1/ cos?] follows 
because F X bb u depends only on the combination of variables (z — C\t) and any 
derivatives of &xBB n will keep the (z — c\t) dependence structure. 

Now, the Poynting vector P X bb u and the energy density Uxbb„ for Fxbb u are 
obtained by considering the real parts of ExBB n and Bxbb„- We have 

(PxbbJ p = -Re{(E XBBn ) e }Re{(B XBBn ) z }; (3.43a) 
(PxBB n )e = Re{(ExBB n ) P }Re{(BxBB n ) z }; (3.43b) 
(PxbbJz = cos i] \\Re{(E XBBn ) P }\ 2 + \Re{(E XBBn )e}\ 2 ] 5 (3.43c) 



v>XBB n = (l + cos 2 r/) [|i2e{fe BBn ) p }| 2 + |i?e{(S XBB J e }| 2 J + \Re{(B XBBn ) z }\ 2 ■ 

(3.44) 

The total energy of F XB B n is then 

/7T /* + 00 /*0O 

d0 / dz / pdpu XB B n (3.45) 
-7T J —CO JO 

Since as z — > 00, E X BB n decreases as l/\z — t cost/ I 1 / 2 , what occurs for the X- 
branches of F X bb„, £xbb„ may not be finite. Nevertheless, as in the case of the 
acoustic X-waves discussed in section 2, we are quite sure that a FAAFxbl„ can 
be launched over a large distance. Obviously in this case the total energy of the 
FAAF XBLn is finite. 

We now restrict our attention to F X bb - In this case from eq.(3.40) and eqs.(3.43) 
we see that (E XBBo ) P = (B XBBo )e = (PxBB )e = 0. In Fig. 4W we see the 
amplitudes of Re{$ X BB } [4(1)], Re{(E XBBo )e} [4(2)], Re{(B XBBo ) p } [4(3)] and 
Re{(B XB B )z} [4(4)]. Fig. 5 shows respectively (P XB b ) p [5(1)], (Pxbb () ) 2 [5(2)] 
and uxbbq [5(3)]. The size of each panel in Figures 4 and 5 is 4m (p-direction) x 

(*)Figures 4, 5 and 6 were reprinted with permission from' 5 !. 



18 



2mm (z-direction) and the maxima and minima of the images in Figures 4 and 5 
(before scaling) are shown in Table 1, in MKSA units^**- 1 . 





Re{$XBB } 


Re{(E XBBo ) e } 


Re{(B XBBo ) p } 


Re{(B XBBo ) z } 


max 


1.0 


9.5 x 10 6 


2.5 x 10 4 


6.1 


min 


0.0 


-9.5 x 10 6 


-2.5 x 10 4 


-1.5 





(Pxbb )p 


(PxBB )z 


Uxbbq 


max 


2.4 x 10 7 


2.4 x 10 11 


1.6 x 10 3 


mm 


-2.4 x 10 7 


0.0 


0.0 



Table 1: Maxima and Minima of the zeroth-order nondiffracting 
electromagnetic X waves (units: MKSA). 

Fig. 6 shows the beam plots of F X bb in Fig. 4 along one of the X-branches 
(from left to right). Fig. 6(1) represents the beam plots of Re{& XBBo } (full line), 

— * — * — * 

Re{(E XBBo ) e } (dotted line), Re{(B XBBo ) p } (dashed line) and Re{(B XBB(] ) z } (long 
dashed line). Fig. 6(2) represents the beam plots of (P XBBo ) p (full line), (P XBB() ) Z 
(dotted line) and u XBBo (dashed line). 

3.4 Finite Aperture Approximation to F XBBo and F XBLo 

From eqs.(3.40), (3.43) and (3.44) we see that E XBBo , B XBB() , P XBBo and u XBBo 
are related to the scalar field & XBBo . It follows that the depth of the field' 5 ' (or non 
diffracting distance — see section 2) of the FAAF XBBo and of the FAAFxbl , which 
of course are to be produced by a finite aperture radiator, are equal and given by 

Z max = £>/2cot77, (3.46) 

where D is the diameter of the radiator and i] is the axicon angle. It can be proved 
also^ that for $xbl ( an d more generally for $xbl„), that Z ma _ x is independent of 
the central frequency of the spectrum B{k) in eq.(3.1). Then if we want, e.g., that 
F XBBo or F XB l travel 115 km with a 20 m diameter radiator, we need rj = 0.005°. 

Figure 7 shows the envelope of Re{FAAQ XBBo } obtained with the finite aperture 
approximation (FAA) given by eq.(A.28), with D = 20 m, a = 0.05 mm and 
rj = 0.005°, for distances z = 10 km [6(1)] and z = 100 km [6(2)], respectively, 

(""^Reprinted with permission from Table I of' 5 '. 
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from the radiator which is located at the plane z — 0. Figures 7(3) and 7(4) show 
the envelope of i?e{FAA$ xbl } f° r the same distances and the same parameters 
(D, a and 77) where B{k) is the following Blackman window function, peaked at the 
frequency f = 700 GHz with a 6 dB bandwidth about 576 GHz: 

m = I a [0.42- 0.5 cos f + 0.08 cos f ], 0<k<2k o] ^ 
\ otherwise; 

where k = 2irf /c (c = 300, 000km/ s). From eq.(3.46) it follows that for the above 
choice of D, a and r\ 

Z max = 115 km (3.48) 

Figs. 8(1) and 8(2) show the lateral beam plots and Figs. 8(3) and 8(4) show the 
axial beam plots respectively for Re{FAA& X BB } an d for Re { FA A <&xbl } used to 
calculate F X bb an d F X bl - The full and dotted lines represent X-waves at distances 
z — 10 km and z = 100 km. Fig. 9 shows the peak values of Re{FAAQ X BB } 
(full line) and Re{FAA<&xBL } (dotted line) along the z-axis from z = 3.45 km to 
z = 230 km. The dashed line represents the result of the exact $xbb solution. The 
6 dB lateral and axial beam widths of <&xbb , which can be measured in Fig 7(1) 
and 7(2), are about 1.96 m and 0.17 mm respectively, and those of the FAA$xbl ( , 
are about 2.5 m and 0.48 mm as can be measured from 7(3) and 7(4). For &xbb 
we can calculate^ 43,26 ! the theoretical values of the 6 dB lateral (BWl) and axial 
(BWa) beam widths, which are given by 

M 1 = ^; BW A = p^. (3.49) 
I sin 77 1 I cos 77I 

With the values of D, ao and 77 given above, we have BWl = 1-98 m and BWa = 
0.17 mm. These are to be compared with the values of these quantities for the 

FAA&xBLo ■ 

We remark also that eq.(3.46) says that Z max does not depend on a . Then we 
can choose an arbitrarily small a to increase the localization (reduced BWl and 
BWa) of the X-wave without altering Z max . Smaller ao requires that the FAA<E» xbl 
be transmitted with broader bandwidth. The depths of field of §xbb and of $xbl 
that we can measure in Fig. 9 are approximately 109 km and 110 km, very close to 
the value given by eq.(3.46) which is 115 km. 

We conclude this section with the following observations. 

(i) In general both subluminal and superluminal UPWs solutions of ME have non 
null field invariants and are not transverse waves. In particular our solutions 
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have a longitudinal component along the z-axis. This result is important 
because it shows that, contrary to the speculations of Evans' 59 ], we do not 
need an electromagnetic theory with a non zero photon-mass, i.e., with F 
satisfying Proca equation in order to have an electromagnetic wave with a 
longitudinal component. Since Evans presents evidence^ of the existence on 
longitudinal magnetic fields in many different physical situations, we conclude 
that the theoretical and experimental study of subluminal and superluminal 
UPW solutions of ME must be continued. 

(ii) We recall that in microwave and optics, as it is well known, the electromag- 
netic intensity is approximately represented by the magnitude of a scalar field 
solution of the HWE. We already quoted in the introduction that Durnin' 16 ' 
produced an optical J -beam, which as seen from eq.(3.1) is related to $xbb 
(®xbl )- If we take into account this fact together with the results of the 
acoustic experiments described in section 2, we arrive at the conclusion that 
subluminal electromagnetic pulses Jo and also superluminal X-waves can be 
launched with appropriate antennas using present technology. 

(iii) If we take a look at the structure of e.g. the FAA$ xbb [eq.(3.40)] plus 
eq.(A.28) we see that it is a "packet" of wavelets, each one traveling with 
speed c. Nevertheless, the electromagnetic X-wave wave that is an interfer- 
ence pattern is such that its peak travels with speed c/ cost] > 1. (This 
indeed happens in the acoustic experiment with c t— > c s , see section 2). Since 
as discussed above we can project an experiment to launch the peak of the 
FAA$xbb fr° m a point z\ to a point Z2, the question arises: Is the existence 
of superluminal electromagnetic waves in conflict with Einstein's Special Rel- 
ativity? We give our answer to this fundamental issue in section 5, but first 
we discuss in section 4 the speed of propagation of the energy associated with 
a superluminal electromagnetic wave. 
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Figure 4: Real part of field components of the exact solution superluminal electro- 
magnetic X-wave at distance z = ct/ cost] (rj = 0.005°, a = 0.05 mm, n = 0). 
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Figure 5: Poynting flux and energy density of the exact solution superluminal elec- 
tromagnetic X-wave at distance z = ctj cost], (rj = 0.005°, a = 0.05 mm, n = 0). 
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Figure 6: (6.1) Beam plots along the X-branches of F X bb for Re{$xRB } or Hertz 

— * — * — * 

potential, Re{(E X BB )d}, Re{(B XBBo ) p }, and Re{(B XB B a )z}- (6.2) Beam plots for 
(Pxbb )p (full line), (P X bb )z (dotted line) and u XBBo (dashed line). 
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Figure 7: 7(1) and 7(2) show the real part of FAA$ X bb at distances z — 10 km 
and z = 100 km from the radiator located at the plane z = with D = 20 m and 
= 0.005°. 7(3) and 7(4) show the real parts of FAA$ X bl f° r the same distances. 
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Figure 8: Beam plots of scalar X-waves (finite aperture). 
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Figure 9: Peak magnitude of X-waves along the z axis. 
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4. The Velocity of Transport of Energy of the UPWs Solu- 
tions of Maxwell Equations 

Motivated by the fact that the acoustic experiment of section 2 shows that the 
energy of the FAA X-wave travels with speed greater than c s and since we found 
in this paper UPWs solutions of Maxwell equations with speeds < v < oo, the 
following question arises naturally: Which is the velocity of transport of the energy 
of a superluminal UPW (or quasi UPW) solution of ME? 

We can find in many physics textbooks (e.gJ 10 ^) and in scientific papers' 4 ^ the 
following argument. Consider an arbitrary solution of ME in vacuum, dF = 0. Then 
if F = E + \B (see eq.(B.17)) it follows that the Poynting vector and the energy 
density of the field are 

P = ExB, u=-{E 2 + B 2 ). (4.1) 
It is obvious that the following inequality always holds: 

I -PI 

v e = — < 1. (4.2) 

Now, the conservation of energy-momentum reads, in integral form over a finite 
volume V with boundary S = dV 

Eq.(4.3) is interpreted saying that § s dS.P is the field energy flux across the surface 
S = dV, so that P is the flux density — the amount of field energy passing through a 
unit area of the surface in unit time. For plane wave solutions of Maxwell equations, 

v £ = 1 (4.4) 

and this result gives origin to the "dogma" that free electromagnetic fields transport 
energy at speed v e — c — 1. 

However v e < 1 is true even for subluminal and superluminal solutions of ME, 
as the ones discussed in section 3. The same is true for the superluminal modified 
Bessel beam found by Bandt 41 J in 1987. There he claims that since v £ < 1 there is 
no conflict between superluminal solutions of ME and Relativity Theory since what 
Relativity forbids is the propagation of energy with speed greater than c. 



28 



Here we challenge this conclusion. The fact is that as is well known P is not 
uniquely defined. Eq(4.3) continues to hold true if we substitute P i— > P + P' with 
V.P' = 0. But of course we can easily find for subluminal, luminal or superluminal 
solutions of Maxwell equations a P' such that 



\P + P'\ 



u 



> 1. (4.5) 



We come to the conclusion that the question of the transport of energy in superlu- 
minal UPWs solutions of ME is an experimental question. For the acoustic superlu- 
minal X-solution of the HWE (see section 2) the energy around the peak area flows 
together with the wave, i.e., with speed c\ = c s j cost] (although the "canonical" 
formula [eq.(2.10)] predicts that the energy flows with v £ < c s ). Since we can see no 
possibility for the field energy of the superluminal electromagnetic wave to travel 
outside the wave we are confident to state that the velocity of energy transport of 
superluminal electromagnetic waves is superluminal. 

Before ending we give another example to illustrate that eq.(4.2) (as is the 
case of eq.(2.10)) is devoid of physical meaning. Consider a spherical conductor in 
electrostatic equilibrium with uniform superficial charge density (total charge Q) 
and with a dipole magnetic moment. Then, we have 

E = Q^-; B = ^(2cosflr + sin#0) (4.6) 

and 

P = ExB = ^-sinO<p, «=^ + ^(3cos^ + l)). (4.7) 

Thus 

\P\ 2CQrsm9 , n „ , n (A Q . 

7 = ^W^TIj^ 0for ^ ' (48) 

Since the fields are static the conservation law eq.(4.3) continues to hold true, as 
there is no motion of charges and for any closed surface containing the spherical 
conductor we have 

' dS.P = 0. (4.9) 



is 

But nothing is in motion! In view of these results we must investigate whether 
the existence of superluminal UPWs solutions of ME is compatible or not with the 
Principle of Relativity. We analyze this question in detail in the next section. 
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To end this section we recall that in section 2.19 of his book Stratton' 19 ! presents 
a discussion of the Poynting vector and energy transfer which essentially agrees with 
the view presented above. Indeed he finished that section with the words: "By this 
standard there is every reason to retain the Poyinting-Heaviside viewpoint until a 
clash with new experimental evidence shall call for its revision." W 

5. Super luminal Solutions of Maxwell Equations and the 
Principle of Relativity 

In section 3 we showed that it seems possible with present technology to launch 
in free space superluminal electromagnetic waves (SEXWs). We show in the follow- 
ing that the physical existence of SEXWs implies a breakdown of the Principle of 
Relativity (PR). Since this is a fundamental issue, with implications for all branches 
of theoretical physics, we will examine the problem with great care. In section 5.1 
we give a rigorous mathematical definition of the PR and in section 5.2 we present 
the proof of the above statement. 

5.1 Mathematical Formulation of the Principle of Relativity and Its 
Physical Meaning 

In Appendix B we define Minkowski spacetime as the triple (M,g,D), where 
M ~ R 4 , g is a Lorentzian metric and D is the Levi-Civita connection of g. 

Consider now Gm, the group of all diffeomorphisms of M, called the manifold 
mapping group. Let T be a geometrical object defined in A C M. The diffeomor- 
phism h G Gm induces a deforming mapping /i,:Th h*T = T such that: 

(i) If / : M D A -> R, then hJ = fo h' 1 : h(A) -> R 

(ii) If T G secT^ s \A) C sec T(M), where T^ S \A) is the sub-bundle of tensors of 

type (r, s) of the tensor bundle T(M) , then 

{h^T)h e Q l * L °ii ■ ■ ■ i h*uj r , h*Xi, . . . , h^X s ) = T e (u;i, . . . , iv r , X 1? . . . , X s ) 
MX,, G T e A, i = 1, . . . , s, Mujj G T*A, j = 1, . . . , r, Ve G A. 

(iii) If D is the Levi-Civita connection and X, Y G sec TM, then 

(KD h * x KY) he h f = (D x Y) e f Ve G M. (5.1) 
(*)Thanks are due to the referee for calling our attention to this point. 
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If {fn = d/dx^} is a coordinate basis for TA and {6> M = dx^} is the corresponding 
dual basis for T*A and if 

T = T^-tO^ ®...®6^®f^®...®f^ (5.2) 

then 

KT = [T^;;f; o h~ l ]K6^ <g> . . . <g> Kd v ° ® hj m <g> . . . <g> ^/ Mr . (5.3) 

Suppose now that A and h(A) can be covered by the local chart (U, if) of the maximal 
atlas of M, and A C. U, h(A) C [/. Let be the coordinate functions associated 
with (U,rj). The mapping 

x" = x^o h- 1 : h(U) -> iR (5.4) 

defines a coordinate transformation i— ► (x'^) if /i(C/) 3 AU Indeed (x'^) 

are the coordinate functions associated with the local chart (V, tp) where h(U) C V 
and [/ fl V 7^ 0. Now, since it is well known that under the above conditions 
h^djdx^ = d/dx'^ and = dx'^, eqs.(5.3) and (5.4) imply that 

(KT) {x , ll) (he) = T {xlt) (e), (5.5) 

where T/ xAt \(e) means the components of T in the chart (a;' 4 ) at the event e G M, 
i.e., T^ M )(e) = T^-^rr^e)) and where T^ 1 -^ r (x' fl (he)) are the components of 
T = hJT in the basis {h^d/dx^ = d/dx'^}, {h*dx^ = dx'^}, at the point h(e). 
Then eq.(5.6) reads 

T££(x*(he))=T£^(x*(e)), (5.6) 

or using eq.(5.5) 

TZ]t{x'»{e)) = (A- 1 )^ . . . A^T^ix'^e)), (5.7) 

where = dx'^/dx a , etc. 

In appendix B we introduce the concept of inertial reference frames 7 G sec TU, 
U C M by 

#(7,7) = 1 and L>7 = 0. (5.8) 

A general frame Z satisfies g(Z, Z) = 1, with DZ ^ 0. If a = ) G secT*£7, it 
holds . 

(£>a) e = a e ® a e + <7 e + w e + -0 e /i e , e E U CM, (5.9) 
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where a = g(A, ), A = DzZ is the acceleration and where ui e is the rotation tensor, 
<r e is the shear tensor, 9 e is the expansion and h e = g\n e where 

T e M = [Z e ] © [H e \. (5.10) 

H e is the rest space of an instantaneous observer at e, i.e. the pair (e, Z e ). Also 
h e (X,Y) = g e (pX,pY), \/X,Y G T e M and p : T e M -> if e . (For the explicit form 
of uo,a,9, see' 60 '). From eqs.(5.9) and (5.10) we see that an inertial reference frame 
has no acceleration, no rotation, no shear and no expansion. 

We introduce also in Appendix B the concept of a (nacs/J). A (nacs/J) (x^) is 
said to be in the Lorentz gauge if x^ , \x = 0, 1, 2, 3 are the usual Lorentz coordinates 
and I = d/dx° G sec TM. We recall that it is a theorem that putting I = eo = 
d/dx°, there exist three other fields G sec TM such that g(e i) e,) = —1, i = 1,2, 3, 
and = d/dx l . 

Now, let be Lorentz coordinate functions as above. We say that t G Gm is 
a Lorentz mapping if and only if 

a>( e ) = AJJa^(e), (5.11) 

where A^ G is a Lorentz transformation. For abuse of notation we denote the 
subset {£} of Gm such that eq.(5.12) holds true also by C\ C Gm- 

When are Lorentz coordinate functions, are also Lorentz coordinate 
functions. In this case we denote 

e M = d/dx^, e'^ = d/dx'^, lfl = dx^, ^ = dx'^ ; (5.12) 

when I G C\ C Gm we say that £*T is the Lorentz deformed version of T. 
Let h G G M - If for a geometrical object T we have 

/i*T = T, (5.13) 

then h is said to be a symmetry of T and the set of all {h G Gm} such that eq.(5.13) 
holds is said to be the symmetry group of T. We can immediately verify that for 
£ G C\ C Gm 

Lg = g, LD = D, (5.14) 

i.e., the special restricted orthochronous Lorentz group C\ is a symmetry group of 
g and D. 

I n [62] we ma intain that a physical theory r is characterized by: 
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(i) the theory of a certain "species of structure" in the sense of Boubarki' 63 ' ; 

(ii) its physical interpretation; 

(iii) its present meaning and present applications. 

We recall that in the mathematical exposition of a given physical theory r, the 
postulates or basic axioms are presented as definitions. Such definitions mean that 
the physical phenomena described by r behave in a certain way. Then, the definitions 
require more motivation than the pure mathematical definitions. We call coordina- 
tive definitions the physical definitions, a term introduced by ReichenbachI 64 ! . It is 
necessary also to make clear that completely convincing and genuine motivations for 
the coordinative definitions cannot be given, since they refer to nature as a whole 
and to the physical theory as a whole. 

The theoretical approach to physics behind (i), (ii) and (iii) above is then to 
admit the mathematical concepts of the "species of structure" defining r as prim- 
itives, and define coordinatively the observation entities from them. Reichenbach 
assumes that "physical knowledge is characterized by the fact that concepts are not 
only defined by other concepts, but are also coordinated to real objects". However, 
in our approach, each physical theory, when characterized as a species of structure, 
contains some implicit geometric objects, like some of the reference frame fields de- 
fined above, that cannot in general be coordinated to real objects. Indeed it would 
be an absurd to suppose that all the infinity of IRF that exist in M must have a 
material support. 

We define a spacetime theory as a theory of a species of structure such that, if 
Mod r is the class of models of r, then each T e Mod r contains a substructure 
called spacetime (ST). More precisely, we have 

T = (ST,T 1 ...T m }, (5.15) 

where ST can be a very general structure ' 6 ^. For what follows we suppose that 
ST = M. = (M,g,D), i.e. that ST is Minkowski spacetime. The Tj, % = 1, . . . ,m 
are (explicit) geometrical objects defined in U C M characterizing the physical fields 
and particle trajectories that cannot be geometrized in T. Here, to be geometrizable 
means to be a metric field or a connection on M or objects derived from these 
concepts as, e.g., the Riemann tensor or the torsion tensor. 

The reference frame fields will be called the implicit geometrical objects of r, 
since they are mathematical objects that do not necessarily correspond to properties 
of a physical system described by r. 
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Now, with the Clifford bundle formalism we can formulate in C£(M) all modern 
physical theories (see Appendix B) including Einstein's gravitational theory'^. We 
introduce now the Lorentz-Maxwell electrodynamics (LME) in C£(M) as a theory 
of a species of structure. We say that LME has as model 

^lme = (M,g,D,F, J, {y9j,m;,e;}), (5.16) 

where (M,g,D) is Minkowski spacetime, {<pi,mi,ei}, i = 1,2, ...,N is the set of 
all charged particles, m 8 and being the masses and charges of the particles and 
tpi : M D / — > M being the world lines of the particles characterized by the fact that 
if (p^ G sec TM is the velocity vector, then (pi = g(tpi*, ) G secA 1 (M) C sec C£(M) 
and (pi.(fi — 1. F G secA 2 (M) C secC£(M) is the electromagnetic field and J G 
sec A X (M) C sec C£(M) is the current density. The proper axioms of the theory are 

OF = J 

rriiD^ipi = eiCpi ■ F 

From a mathematical point of view it is a trivial result that tlme has the 
following property: If h G Gm and if eqs.(5.16) have a solution (F, J, (tp i7 e^)) in 
U C. M then {h*F,h*J, (h*(pi,m,i,ei)) is also a solution of eqs.(5.16) in h(U). Since 
the result is true for any h G Gm it is true for £ G C Gm, i-e., for any Lorentz 
mapping. 

We must now make it clear that (F, J, {{pi, rrii, e^}) which is a solution of eq.(5.16) 
in U can be obtained only by imposing mathematical boundary conditions which we 
denote by BU . The solution will be realizable in nature if and only if the mathe- 
matical boundary conditions can be physically realizable. This is indeed a nontrivial 
point' 62 ] for in particular it says to us that even if {h*F, h*J, {h*(fi, e^}) can be a 
solution of eqs.(5.16) with mathematical boundary conditions Bh(U), it may hap- 
pen that Bh(U) cannot be physically realizable in nature. The following statement, 
denoted PR±, is usually presented' 62 ^ as the Principle of (Special) Relativity in active 
form: 

Pi?i:Let £ G C + C Gm- If for a physical theory r and T G Modr, T = 
(M, g, D, Ti, . . . , T m ) is a possible physical phenomenon, then £*Y = (M, g, D, 
Z*Ti, . . . , /*T m ) is also a possible physical phenomenon. 

It is clear that hidden in PRi is the assumption that the boundary conditions 
that determine £*Y are physically realizable. Before we continue we introduce the 
statement denoted PR2, known as the Principle of (Special) Relativity in passive 

f orm [62] 
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(5.17) 



Pi?2:"All inertial reference frames are physically equivalent or indistinguishable". 



We now give a precise mathematical meaning to the above statement. 

Let r be a spacetime theory and let ST = (M, g, D) be a substructure of Mod r 
representing spacetime. Let I G sec TU and I' G sec TV , U, V C M, be two inertial 
reference frames. Let {U,rj) and (V,(p) be two Lorentz charts of the maximal atlas 
of M that are naturally adapted respectively to I and V . If and (x M ) are the 
coordinate functions associated with (U,r]) and (V, </?), we have I = d/dx°,I' = 
d/dx °. 

Definition: Two inertial reference frames / and I' as above are said to be phys- 
ically equivalent according to r if and only if the following conditions are satisfied: 

(i) G M D C\ 3 £ : U -> £{U) C V, a> = ac" o r 1 I' = £J 

When T G Modr, T = (M, g>, -D, Ti, . . . T m ), is such that g and D are defined 
over all M and T< G secC£(ff) C sec C£(M), calling o = (g, D, T h ... T m ), o solves 
a set of differential equations in rj(U) C M 4 with a given set of boundary conditions 
denoted b°^\ which we write as 

D? x » ) (o {x » ) ) e = ; ■ eeU (5.18) 

and we must have: 

(ii) If T G Mod t £*T G Mod r, then necessarily 

4T = (M, g, D, £*Ti, . . . £*T m ) (5.19) 
is defined in £(U) C V and calling £*o = {g, D, £*Ti, . . . , £*T m } we must have 

D " x ',>)(e*°(x'*))\te = ; be * o{x ' n £e e *(u) £ v - ( 5 - 2 °) 



In eqs.( |5.18| ) and ( |5.2U| ) D?^ and D? x <^ mean a = 1,2, ... , m sets of differential 



equations in iR 4 . The system of differential equations (5.19) must have the same 

functional form as the system of differential equations (5.17) and b e * o( ~ x M ^ must be 
relative to (x M ) the same as b ^ is relative to (x^) and if b°( x>i < is physically realiz- 
able then b e *°( x M ^ must also be physically realizable. We say under these conditions 
that I ~ I' and that £*o is the Lorentz deformed version of the phenomena described 
by o. 

Since in the above definition £ if T = (M, g,D,£*Ti, ... ,£*T m ), it follows that 
when / ~ I', then £*g = g,£*D = D (as we already know) and this means that the 
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spacetime structure does not give a preferred status to I or /' according to r. 

5.2 Proof that the Existence of SEXWs Implies a Breakdown of Piix 
and PR 2 

We are now able to prove the statement presented at the beginning of this sec- 
tion, that the existence of SEXWs implies a breakdown of the Principle of Relativity 
in both its active (PRi) and passive (PR2) versions. 

Let i e L\ C G m and let F, F E secA 2 (M) C sec C£(M), F = £*F. Let F = 
£*F = RFR~\ where F e = (l/2)F fll/ (x s (£- 1 e))^Y and where R £ sec Spin + (1, 3) C 
secC£(M) is a Lorentz mapping, such that 7'^ = R'-f^R^ 1 = A^7°, e C\ and let 
(x M ) and (x M ) be Lorentz coordinate functions as before such that 7^ = <ix M ,7 M = 
dx M and x M = x M o £~ l . We write 

F e = ^ u (x 6 (e))rr; (5.22a) 

F e = ^(^(e)hV'; (5.22b) 

F e = ±F„ v (x 5 (e))rr; (5.23a) 

F e = ^(x' 5 (e))7 V 7 V . (5.23b) 
From (5.22a) and (5.22b) we get that 

FLp(x' S (e)) = {h-- 1 Y a {h--%F, v {x\e)). (5.24) 
From (5.22a) and (5.23b) we also get 

F aP {x\e)) = K^F, u {x\r l e)) (5.25) 

Now, suppose that F is a superluminal solution of Maxwell equation, in par- 
ticular a SEXW as discussed in section 3. Suppose that F has been produced 
in the inertial frame / with (x^) as (nacs/J), with the physical device described 
in section 3. F is generated in the plane z = and is traveling with speed 
c\ — \j cos r] in the negative z-direction. It will then travel to the future in 
spacetime, according to the observers in /. Now, there exists £ G C\ such that 
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£*F = F = RFR~ l will be a solution of Maxwell equations and such that if the 
velocity 1-form of F is vp = (c\ — \)~ 1 I 2 {1, 0, 0, — C\), then the velocity 1-form of 
F is v p = (c'i — l) _1 / 2 (— 1, 0, 0, — c^), with c[ > 1, i.e., v-p is pointing to the past. 
As its is well known F carries negative energy according to the observers in the / 
frame. 

We then arrive at the conclusion that to assume the validity of PR\ is to assume 
the physical possibility of sending to the past waves carrying negative energy. This 
seems to the authors an impossible task, and the reason is that there do no exist 
physically realizable boundary conditions that would allow the observers in I to 
launch F in spacetime and such that it traveled to its own past. 

We now show that there is also a breakdown of PR2, i-e., that it is not true that 
all inertial frames are physically equivalent. Suppose we have two inertial frames / 
and V as above, i.e., I = d/dx°, V = d/dx °. 

Suppose that F is a SEXW which can be launched in I with velocity 1-form as 
above and suppose F is a SEXW built in I' at the plane z' = and with velocity 
1-form relative to (x M ) given by v-p = v and 

^=( T =L=,0,0,- T |L=) (5.26) 

If F and F are related as above we see (See Fig. 10) that F, which has positive 
energy and is traveling to the future according to I', can be sent to the past of the 
observers at rest in the / frame. Obviously this is impossible and we conclude that 
F is not a physically realizable phenomenon in nature. It cannot be realized in V 
but F can be realized in /. It follows that PR2 does not hold. 

If the elements of the set of inertial reference frames are not equivalent then there 
must exist a fundamental reference frame. Let I G sec TM be that fundamental 
frame. If V is moving with speed V relative to /, i.e., 

1 T / 

/' = tt ~ ; d/dz , (5.27) 

then, if observers in V are equipped with a generator of SEXWs and if they prepare 
their apparatus in order to send SEXWs with different velocity 1-forms in all pos- 
sible directions in spacetime, they will find a particular velocity 1-form in a given 
spacetime direction in which the device stops working. A simple calculation yields 
then, for the observes in the value of V ! 

In ' 65 1 Recami argued that the Principle of Relativity continues to hold true even 
though superluminal phenomena exist in nature. In this theory of tachyons there 
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exists, of course, a situation completely analogous to the one described above (called 
the Tolman-Regge paradox), and according to Recami's view PR2 is valid because 
V must interpret F a being an anti-SEXW carrying positive energy and going into 
the future according to him. In his theory of tachyons Recami was able to show that 
the dynamics of tachyons implies that no detector at rest in / can detect a tachyon 
(the same would be valid for a SEXW like F) sent by V with velocity 1-form given 
by eq.(4.26). Thus he claimed that PR2 is true. At first sight the argument seems 
good, but it is at least incomplete. Indeed, a detector in I does not need to be at 
rest in J. We can imagine a detector in periodic motion in / which could absorb 
the F wave generated by V if this was indeed possible. It is enough for the detector 
to have relative to I the speed V of the /' frame in the appropriate direction at the 
moment of absorption. This simple argument shows that there is no salvation for 
PR2 (and for PR\) if superluminal phenomena exist in nature. 

The attentive reader at this point probably has the following question in his/her 
mind: How could the authors start with Minkowski spacetime, with equations car- 
rying the Lorentz symmetry and yet arrive at the conclusion that PR\ and PR2 do 
not hold? 

The reason is that the Lorentzian structure of (M,g,D) can be seen to exist 
directly from the Newtonian spacetime structure as proved in t 66 l. In that paper 
Rodrigues and collaborators show that even if C\ is not a symmetry group of New- 
tonian dynamics it is a symmetry group of the only possible coherent formulation of 
Lorentz-Maxwell electrodynamic theory compatible with experimental results that 
is possible to formulate in the Newtonian spacetime^. 

We finish calling to the reader's attention that there are some experiments 
reported in the literature which suggest also a breakdown of PR 2 for the roto- 
translational motion of solid bodies. A discussion and references can be found int 67 l 

6. Conclusions 

In this paper we presented a unified theory showing that the homogeneous wave 
equation, the Klein-Gordon equation, Maxwell equations and the Dirac and Weyl 
equations have solutions with the form of undistorted progressive waves (UPWs) of 
arbitrary speeds < v < 00. 

We present also the results of an experiment which confirms that finite aperture 
approximations to a Bessel pulse and to an X-wave in water move as predicted by 

(*) We recall that Maxwell equations have, as is well known, many symmetry groups besides 
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Figure 10: F cannot be launched by V . 



our theory, i.e., the Bessel pulse moves with speed less than c s and the X-wave 
moves with speed greater than c s , c s being the sound velocity in water. 

We exhibit also some subluminal and superluminal solutions of Maxwell equa- 
tions. We showed that subluminal solutions can in principle be used to model 
purely electromagnetic particles. A detailed discussion is given about the superlu- 
minal electromagnetic X-wave solution of Maxwell equations and we showed that it 
can in principle be launched with available technology. Here a point must be clear, 
the X-waves, both acoustic and electromagnetic, are signals in the sense defined by 
Nimtzt 74 '. It is a widespread misunderstanding that signals must have a front. A 
front can be defined only mathematically because it implies an infinite frequency 
spectrum. Every real signal does not have a well defined front. 

The existence of superluminal electromagnetic waves implies in the breakdown 
of the Principle of Relativity. W We observe that besides its fundamental theoretical 
implications, the practical implications of the existence of UPWs solutions of the 
main field equations of theoretical physics (and their finite aperture realizations) are 

Wit is important to recall that there exists the possibility of propagation of superluminal sig- 
nals inside the hadronic matter. In this case the ingenious construction of Santilli's isominkowskian 
spaces (seet 68_73 l) is useful. 
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very important. This practical importance ranges from applications in ultrasound 
medical imaging to the project of electromagnetic bullets and new communication 
devices^ 33 '. Also we would like to conjecture that the existence of subluminal and 
superluminal solutions of the Weyl equation may be important to solve some of the 
mysteries associated with neutrinos. Indeed, if neutrinos can be produced in sublu- 
minal or superluminal modes — see^ 75,76 ' for some experimental evidence concerning 
superluminal neutrinos — they can eventually escape detection on earth after leaving 
the sun. Moreover, for neutrinos in a subluminal or superluminal mode it would be 
possible to define a kind of "effective mass" . Recently some cosmological evidences 
that neutrinos have a non- vanishing mass have been discussed by e.g. Primack et 
al' 77 ]. One such "effective mass" could be responsible for those cosmological evi- 
dences, and in such a way that we can still have a left-handed neutrino since it 
would satisfy the Weyl equation. We discuss more this issue in another publication. 
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Appendix A. Solutions of the (Scalar) Homogeneous Wave 
Equation and Their Finite Aperture Realizations 

In this appendix we first recall briefly some well known results concerning the 
fundamental (Green's functions) and the general solutions of the (scalar) homoge- 
neous wave equation (HWE) and the theory of their finite aperture approximation 
(FAA). FAA is based on the Rayleigh-Sommerf eld formulation of diffraction (RSFD) 
by a plane screen. We show that under certain conditions the RSFD is useful for 
designing physical devices to launch waves that travel with the characteristic veloc- 
ity in a homogeneous medium (i.e., the speed c that appears in the wave equation). 
More important, RSFD is also useful for projecting physical devices to launch some 
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of the subluminal and superluminal solutions of the HWE (i.e., waves that propa- 
gate in an homogeneous medium with speeds respectively less and greater than c) 
that we present in this appendix. We use units such that c = 1 and H — 1, where 
c is the so called velocity of light in vacuum and Ti is Planck's constant divided by 2n. 

Al. Green's Functions and the General Solution of the (Scalar) HWE 

Let $ in what follows be a complex function in Minkowski spacetime M: 

$:I9ih G <F . (A.l) 

The inhomogeneous wave equation for $ is 

D$ = (Jj£ - V2 ) $ = ' ( A - 2 ) 

where p is a complex function in Minkowski spacetime. We define a two-point 
Green's function for the wave equation (A.2) as a solution of 

DG(x - x) = 4ir5(x - x') . (A.3) 

As it is well known, the fundamental solutions of (A.3) are: 

Retarded Green's function: Gr(x — x') = 2H(x — x')5[(x — x') 2 }; (A. 4a) 
Advanced Green's function: G A (x - x') = 2H[-(x - x')]S[(x - x') 2 ]; (A.4b) 

where (x — x') 2 = (x° — x' ) 2 — (x — x') 2 , H(x) = H(x°) is the step function and 
x o = tjX 'o = 1 J. 

We can rewrite eqs.(A.4) as (R — \x — x'\): 

G R {x° - x'°; x-x') = ]-5(x° - x'° - R) ; (A.4c) 

R 

G A (x° - x'°; x-x') = ]-5{x Q - x'° - R) . (A.4d) 

R 

We define the Schwinger function by 

G S = G R -G A = 2e{x)5{x 2 ); e{x) = H{x) - H{-x) . (A.5) 
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It has the properties 



nG s = 0; G s (x) = -G s (-x); G s (x) = if x 2 < ; 



G s (0,x) = 0; 



dG f 



dx i 



0; 



dG H 



x u =0 



S(x). 



(A.6a) 
(A.6b) 



For the reader who is familiar with the material presented in Appendix B, we 
observe that these equations can be rewritten in a very elegant way in C£c(M). (If 
you haven't read Appendix B, go to eq.(A.8').) We have 



/ *dG s {x -y) = - j dG s {x - y)7 5 = 1, if y E a, 

J a J a 



(A.7) 



where a is any spacelike surface. Then if / e sec (& ® /\°(M) C secC£c(M) is any 
function defined on a spacelike surface a, we can write 

J[*dG 8 {x - y)]f(x) = -J dG s (x - y)f(xh 5 = f(y) . (A.8) 
Eqs.(A.7) and (A.8) appear written in textbooks on field theory as 

/ d,G s (x - y)da^(x) = 1 ; / f(x)d,G s (x - y)da^(x) = f(y) . (A.8') 



We now express the general solution of eq.(A.2), including the initial conditions, in 
a bounded constant time spacelike hypersurface a characterized by 7 1 A 7 2 A 7 s in 
terms of Gr. We write the solution in the standard vector notation. Let the constant 
time hypersurface a be the volume V C 1R 3 and dV = S its boundary. We have, 



$(t, x) = J* + dt' J jj dv'G R (t -t',x- x')p(t', x') 



+ 



1 

1 r 



9$ 



d 



G R \ t , =0 —{t',x!)\ t , =Q - ®(t',x%, =0 —G R \ t , =0 
~t Io + dt ' Sis rf ^'-( G '«g rad,$ - ®%™dG R ), 



(A.9) 



where grad' means that the gradient operator acts on x', and where t + means that 
the integral over t' must end on t' = t+e in order to avoid ending the integral exactly 
at the peak of the 5-function. The first term in eq.(A.9) represents the effects of 
the sources, the second term represents the effects of the initial conditions (Cauchy 
problem) and the third term represents the effects of the boundary conditions on 
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the space boundaries dV = S.This term is essential for the theory of diffraction and 
in particular for the RSFD. 

d 

Cauchy problem: Suppose that $(0, x) and Tr&(t, x)\ t=0 are known at every point 

in space, and assume that there are no sources present, i.e., p = 0. Then the solution 
of the HWE becomes 



G R \ t , =0 -<f>(t',x%, =0 - -GU, =0 $(0,f ) 
at at 



(A.10) 



The integration extends over all space and we explicitly assume that the third term 
in eq.(A.9) vanishes at infinity. 

We can give an intrinsic formulation of eq.(A.lO). Let x G a, where a is a 
spacelike surface without boundary. Then the solution of the HWE can be written 



(A.ll) 



$(ar) = — / {G s (x - x')[*d$(x')} - [*dG s (x - x')]<S>(x')} 

4:71 Ja 

= ^ Jj^(x)[G s (x - x')d^(x') - d,G s (x - x')$(x')} 

where Gs is the Schwinger function [see eqs.(A.7, A. 8)]. <&(x) given by eq.(A.ll) 
corresponds to "causal propagation" in the usual Einstein sense, i.e., <J>(x) is in- 
fluenced only by points of a which lie in the backward (forward) light cone of x' , 
depending on whether x is "later" ("earlier") than a. 

A2. Huygen's Principle; the Kirchhoff and Rayleigh-Sommerfeld Formu- 
lations of Diffraction by a Plane Screen' 79 ' 



Huygen's principle is essential for understanding Kirchhoff 's formulation and the 
Rayleigh-Sommerfeld formulation (RSF) of diffraction by a plane screen. Consider 
again the general solution [eq.(A.9)j of the HWE which is non-null in the surface 

d 

S = dV and suppose also that $(0, x) and ir&(t, x)\ t= o are null for all x G V. Then 

CJ v 



eq.(A.9) gives 



^M(t'^) + ^(tW)-^Ht', 



i? 3 



R 2 dt> 



' x') 



t'=t-R 

(A.12) 
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From eq.(A.12) we see that if S is along a wavefront and the rest of it is at infinity 
or where <£> is zero, we can say that the field value <3> at (t, x) is caused by the field 
$ in the wave front at time (t — R) earlier. This is Huygen's principle. 

Kirchhoff's theory: Now, consider a screen with a hole like in Fig.ll. 



Figure 11: Diffraction from a finite aperture. 
Suppose that we have an exact solution of the HWE that can be written as 

$(t,£) = F{x)e iuj \ (A. 13) 

where we define also 

u = k (A. 14) 

and k is not necessarily the propagation vector (see bellow). We want to find the 
field &txeV, with dV = S 1 + S 2 (Fig.ll), with p = Vf e V. Kirchhoff proposed 
to use eq.(A.12) to give an approximate solution for the problem. Under the so 
called Sommerfeld radiation condition, 

lim r I ^ - ikF } = 0, (A. 15) 

where r = \f\ = x — x', x' being a point of S 2 , the integral in eq.(A.12) is null over 
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S 2 . Then, we get 



— ikR 

G K = — — — , R = \x — x*\, x 1 £ Si . (A.17) 
R 

Now, the "source" is opaque, except for the aperture which is denoted by E in 
Fig.ll. It is reasonable to suppose that the major contribution to the integral arises 
from points of Si in the aperture S C Si. Kirchhoff then proposed the conditions: 

(i) Across E, the fields F and dF/dn are exactly the same as they would be in 
the absence of sources. 

(ii) Over the portion of Si that lies in the geometrical shadow of the screen the 
field F and dF/dn are null. 

Conditions (i) plus (ii) are called Kirchhoff boundary conditions, and we end 
with 

F *®=JL*{i&*- F &*)- (A18) 

where Fk(x) is the Kirchhoff approximation to the problem. As is well known, Fk 
gives results that agree very well with experiments, if the dimensions of the aperture 
are large compared with the wave length. Nevertheless, Kirchhoff's solution is in- 
consistent, since under the hypothesis given by eq.(A.13), F(x) becomes a solution 
of the Helmholtz equation 

V 2 F + cj 2 F = 0, (A.19) 

and as is well known it is illicit for this equation to impose simultaneously arbitrary 
boundary conditions for both F and dF/dn. 

A further shortcoming of Fx is that it fails to reproduce the assumed bound- 
ary conditions when x e E C Si. To avoid such inconsistencies Sommerfeld pro- 
posed to eliminate the necessity of imposing boundary conditions on both F and 
dF/ dn simultaneously. This gives the so called Rayleigh- Sommerfeld formulation of 
diffraction by a plane screen (RSFD). RSFD is obtained as follows. Consider again 
a solution of eq.(A.18) under Sommerfeld radiation condition [eq.(A.15)j 

^-jis^-^)^ <A2o) 
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where now Grs is a Green function for eq.(A.19) different from Gk- Grs must 
provide an exact solution of eq.(A.19) but we want in addition that Grs or dGRs/dn 
vanish over the entire surface Si, since as we already said we cannot impose the 
values of F and dF/dn simultaneously. 

A solution for this problem is to take Grs as a three-point function, i.e., as a 
solution of 



(V 2 + oj 2 )Gr S (x , f , x") = An5{x - f ) - 47T(5(f - f ')■ 



We get 



G RS (x, x', x") j 

| 3v 3/ I . I 3/ 3/ I . 



(A.21) 

(A.22) 
(A.23) 



where x G Si and x' = —x" are mirror image points relative to Si. This solution 



gives Grs 



and dG RS /dn 



Si 



7^0. 



Si 



7^ and dG^ s /dn 



Another solution for our problem such that G~^ s 
realized for Gr~ s satisfying 

(V 2 + uj 2 )G + RS {x , x' , x") = Att5(x - x') + Aix5(x - x"). 

Then 



is 



Si 



Grs(%i x ,x') 



^ikR gikR' 
+ 



(A.24) 



(A.25) 



R R' 
with R and R' as in eq.(A.23). 

We now use Gr S m eq.(A.25) and take Si as being the z = plane. In this case 
n = —k, k being the versor of the z direction, R = x — x', R.n = z' — z cos(n, R) = 
(z' — z)/R and we get 



f ^ = ~tJ L iS ' F(x '> y '> a) 



(A.26) 



A3. Finite Aperture Approximation for Waves Satisfying x) = F(x)e 



-iu:t 



The finite aperture approximation to eq.(A.26) consists in integrating only over 
£ C Si, i.e., we suppose F(x) = Vrr G (Si\S). Taking into account that 

k = 2tt/A, u = k, (A.27) 
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we get 



In section A4 we show some subluminal and superluminal solutions of the HWE 
and then discuss for which solutions the FAA is valid. We show that there are indeed 
subluminal and superluminal solutions of the HWE for which (A. 28) can be used. 
Even more important, we describe in section 2 the results of recent experiments, 
conducted by us, that confirm the predictions of the theory for acoustic waves in 
water. 



A4. Subluminal and Superluminal Solutions of the HWE 



Consider the HWE (c = 1) 

$ - V 2 $ = . (A.2') 

at 1 

We now present some subluminal and superluminal solutions of eq.(A.2') J 80 ' 

Subluminal and Superluminal Spherical Bessel Beams. To introduce these beams we 
define the variables 

U = [x 2 + y 2 + ll(z-v < tf}^- (A.29a) 

7< = ^— ; o£-**=n 2 <; v < = ^; (A.29b) 
J I — v< " ft < 

i> = [-x 2 -y 2 + i 2 > {z-v > tf] 1 ' 2 - (A.29c) 

7> = . 1 ; ul - k> = -ft> ; v> = du > /dk > . (A.29d) 
\l v l - 1 

We can now easily verify that the functions and <3>^™ below are respectively 
subluminal and superluminal solutions of the HWE (see example 3 below for how 
to obtain these solutions). We have 

¥™(t,x) = C iJe (n p Q P^cosOy^e 1 ^-^) (A.30) 

where the index p =<, >, Ce are constants, je are the spherical Bessel functions, 
P^ are the Legendre functions and (r, 9, ip) are the usual spherical coordinates. 
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§ e ™ has phase velocity (w < /k < ) < 1 [(w>/k>) > 1] and the modulation 

function je(Q<£<) \je(^>C>)] moves with group velocity t>< [i>>], where < v < < 1 
[1 < v > < oo]. Both and are undistorted progressive waves (UPWs). This 
term has been introduced by Courant and HilbertW; however they didn't suspect of 
UPWs moving with speeds greater than c = 1. For use in the main text we write 
the explicit form of $<° and $° > °, which we denote simply by $< and $>: 

%(t,x)=C Sm{ ^ p) e l ^ t ~ k ^-, p=< or > . (A.31) 
When v< = 0, we have $< — > $ , 

$ o(t) f) = c ^£i; e *M r = {x 2 +y 2 + z 2y/2 (A>32) 

r 

When v> = oo, cu> = and $° — > $oo, 

^(t,x) = Cj^e in > z , p = {x 2 + y 2 ) 1 ' 2 . (A.33) 

We observe that if our interpretation of phase and group velocities is correct, 
then there must be a Lorentz frame where <£>< is at rest. It is trivial to verify that 
in the coordinate chart (x M ) which is a 

(nacs/T), where /'=(!- vlY^d/dt + 



(i></ y 1 — v\)d/dz is a Lorentz frame moving with speed v < in the z direction rel- 
ative to / = d/dt, $ p goes in $ (t',x') given by eq.(A.32) with 1 1— > t', a 7 1— > a 7 '. 

Subluminal and Superluminal Bessel Beams. The solutions of the HWE in cylindrical 
coordinates are well known' 19 ]. Here we recall how these solutions are obtained in 
order to present new subluminal and superluminal solutions of the HWE. In what 
follows the cylindrical coordinate functions are denoted by (p, 9, z), p = (x 2 + y 2 ) 1 ^ 2 , 
x = p cos 9, y = psin6>. We write for $: 

$(t,p,9,z) = f 1 (p)f 2 (9)f 3 (t,z). (A.34) 

Inserting (A.34) in (A. 2') gives 

?j? h + % fl + {Bp2 ~ " 2)/l = 0i (A - 35a) 

(^ + ^)A = 0; (A.35b) 
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In these equations B and v are separation constants. Since we want $ to be periodic 
in 9 we choose c = nan integer. For B we consider two cases: 

(i) Subluminal Bessel solution, B = Q< > 

In this case (A. 35a) is a Bessel equation and we have 

$< (t, p, 9, z) = C n J n (P^<y (k<z ~ w<t+n9) , n = 0, 1, 2, . . . , (A.36) 

where C n is a constant, J n is the n-th order Bessel function and 

u%-h% = n 2 < . (A.37) 

I n [43] <j><:^ are ca u e d the nth-order non- diffracting Bessel beams**) . 

Bessel beams are examples of undistorted progressive waves (UPWs). They are 
"subluminal" waves. Indeed, the group velocity for each wave is 

w< = duj < /dk < , < v< < 1 , (A.38) 

but the phase velocity of the wave is (cu < /k < ) > 1. That this interpretation is correct 
follows from the results of the acoustic experiment described in section 2. 

It is convenient for what follows to define the variable 77, called the axicon 
angled, 

k < = k < cosr] , fi< = k< sin 77 , < rj < ir/2 . (A. 39) 

Then 

k < = uj < > (A.40) 
and eq.(A.36) can be rewritten as = with 

$< n = C n J n (k<p sin ^e^^r^M. (A.41) 

In this form the solution is called in' 43 ] the n-th order non- diffracting portion of 
the Axicon Beam. The phase velocity v ph = 1/ cost? is independent of fc<, but, of 
course, it is dependent on k < . We shall show below that waves constructed from the 
$j n beams can be subluminal or superluminal ! 

(ii) Superluminal (Modified) Bessel Solution, B = < 

In this case (A. 35a) is the modified Bessel equation and we denote the solutions by 

p, 0, z) = C n K n {n >P )e^> z ^> t+ne \ n = 0, 1, . . . , (A.42) 

(*)Thc only difference is that fc< is denoted by f3 = \Juj\ — fi< and u> < is denoted by k' = 
ui/c > 0. (We use units where c = 1). 
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where K n are the modified Bessel functions, C n are constants and 

c4 - = -Ml • (A.43) 

We see that $^ are also examples of UPWs, each of which has group velocity 
v > = du > /dk > such that 1 < v > < oo and phase velocity < (u;>//c>) < 1. As in 
the case of the spherical Bessel beam [eq.(A.31)j we see again that our interpretation 
of phase and group velocities is correct. Indeed, for the superluminal (modified) 
Bessel beam there is no Lorentz frame where the wave is stationary. 

The $^- Q beam was discussed by Band' 41 ] in 1988 as an example of superluminal 
motion. Band proposed to launch the beam in the exterior of a cylinder of 
radius r\ on which there is an appropriate superficial charge density. Since _K" (^> r i) 
is non singular, his solution works. In section 3 we discuss some of Band's statements. 

We are now prepared to present some other very interesting solutions of the 
HWE, in particular the so called X-waves' 43 ^ which are superluminal, as proved by 
the acoustic experiments described in section 2. 

Theorem [Lu and GreenleafY 43 ^ : The three functions below are families of exact 
solutions of the HWE [eq.(A.2')] in cylindrical coordinates: 



1 

2^ 



A(4>)f(s)d<p 



roo 

%(s) = / T{k< 
Jo 

* K (s) = f^D( v ) [i- f_^A{<t>)f{s)d4> 



dk< ; (A.44) 
dr] ; (A.45) 



$ L (p, 9,z-t) = $!(p, 6)<& 2 (z - t) ; (A.46) 
s = a (k < ,r])pcos(9 - 0) + b(k <: rj)[z ± c 1 (k <: rj)t] (A.47) 



where 
and 

ci(k<,v) = + 77) /&(£<, r?)] 2 . (A.48) 

In these formulas T(k K ) is any complex function (well behaved) of k < and could 
include the temporal frequency transfer function of a radiator system, A(<f>) is any 
complex function (well behaved) of (f> and represents a weighting function of the inte- 
gration with respect to <p, f(s) is any complex function (well behaved) of s (solution 
of eq.(A.29)), D(rj) is any complex function (well behaved) of rj and represents a 
weighting function of the integration with respect to 77, called the axicon angle (see 
eq.(A.39)), ao(k < ,r)) is any complex function of k K and 77, 6( A;< , 77) is any complex 
function of k < and rj. 
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As in the previous solutions, we take c = 1. Note that k < , 77 and the wave vector 
k < of the f(s) solution of eq.(A.29) are related by eq.(A.39). Also $ 2 (z — t) is any 
complex function of (z — t) and 6) is any solution of the transverse Laplace 
equation, i.e., 



The proof is obtained by direct substitution of $> v , § K and $ L in the HWE. Ob- 
viously, the exact solution $i is an example of a luminal UPW, because if one 
"travels" with the speed c = 1, i.e., with z — t — const., both the lateral and axial 
components 9) and $2(2 — t) will be the same for all time t and distance z. 
When c\{k, 77) in eq.(A.47) is real, (±) represent respectively backward and forward 
propagating waves. 

We recall that & v (s) and <&k{ s ) represent families of UPWs if Ci(k < ,rj) is inde- 
pendent of k < and 77 respectively. These waves travel to infinity at speed c±. & v {s) 
is a generalized function that contains some of the UPWs solutions of the HWE 
derived previously. In particular, if T(k < ) = S(k < — fc < ), k K = to > is a constant 



and if f(s) = e s , ao(k < ,r)) = —ifl < , b{k <) rf) — i/3 — iu/ci, one obtains Durnin's 
UPW beamt 10 ] 



If A((f>) = i n e m<?i we obtain the n-th order UPW Bessel beam $j n given by eq.(A.36). 
$^ n (s) is obtained in the same way with the transformation fc< = k < cosrj; Q < = 
k K sin rj. 

The X-waves. We now present a superluminal UPW wave which, as discussed 
in section 2, is physically realizable in an approximate way (FAA) in the acoustic 
case and can be used to generate Hertz potentials for the electromagnetic field (see 
section 3). We take in eq.(A.44): 




(A.49) 




(A.50) 



T(k < ) = 5(A; < )e- a()fc <; A{<j>) = i n e 
b(k < , rj) = ik cosrj; f(s) = e s . 



vru 



ao(k < ,r]) = —ik < sinr]; 



(A.51) 



We then get 




(A.52) 
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In eq.(A.52) B(k <) is any well behaved complex function of fc< and represents a 
transfer function of practical radiator, k < = to and ao is a constant, and rj is again 
called the axicon angle' 26 '. Eq.(A.52) shows that $^ is represented by a Laplace 
transform of the function B{k < )J n {k < ps\.nri) and an azimuthal phase term e md . 
The name X-waves for the $y comes from the fact that these waves have an X- 
like shape in a plane containing the axis of symmetry of the waves (the z-axis, see 
Fig. 4(1) in section 3). 

The &xBB n waves. This wave is obtained from eq.(A.44) putting B{k < ) = oq. It 
is called the X-wave produced by an infinite aperture and broad bandwidth. We use 
in this case the notation &xBB n - Under these conditions we get 

= ao^sin,)^ (n = o>lj2> .. (A.53) 

where the subscript denotes "broadband" . Also 

M = (pshiT?) 2 + r 2 ; (A.54) 

r = [a — i(z cosr] — t)] (A. 55) 

For n = we get &xbb : 

* >XBB ° = It ■ ,2 r a ° v • (A ' 56) 

J(p smr]) z + [ao — t{z cos 77 — t)\ z 

It is clear that all $j BBn are UPWs which propagate with speed C\ — 1/ cos 77 > 1 
in the z-direction. Our statement is justified for as can be easily seen (as in the 
modified superluminal Bessel beam) there is no Lorentz frame where $^ BBn is at 
rest. Observe that this is the real speed of the wave; phase and group velocity 
concepts are not applicable here. Eq( |A.56p does not give any dispersion relation. 

The $j BB waves cannot be produced in practice as they have infinite energy 
(see section A7), but as we shall show a good approximation for them can be realized 
with finite aperture radiators. 

A5. Construction of $j n and X- Waves with Finite Aperture Radiators 

In section A3 we study the condition under which the Rayleigh-Sommerfeld 
solution to HWE [eq.(A.24)j can be derived. The condition is just that the wave $ 
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must be written as x) = F(x)e~ luJt ; which is true for the Bessel beams In 
section 2 we show that a finite aperture approximation (FAA) to a broad band Bessel 
beam or Bessel pulse denoted FAA§ BBJn or $faj„ [see eq.(2.3)] can be physically 
realized and moves as predicted by the theory. 

At first sight it is not obvious that for the $x„ waves we can use eq.(A.26), but 
actually we can. This happens because we can write, 

1 r+oo ~> 

**„M = teJ-oo**^'® 6 " (A,57) 

$ Xn (u,x) = 2ne me B{uj)J n {ujpsmr])H{u)e'^ ao - izco ^\ n = 0, 1, 2, . .(A.58) 

where H(uS) is the step function and each &x n {^,x) is a solution of the transverse 
Helmholtz equation. Then the Rayleigh-Sommerfeld approximation can be written 
and the FAA can be used. Denoting the FAA to by ^pAx n an d using eq.(A.28) 
we get 

1 p2n i-D/2 ~ _ e ik < R 

® > FAx n (k < ,x) = -J Q d9'j o p>dp'<S> Xn (k<,p', — z 

1 /-27T rD/2 ~ e ik<R 

+ ^J o d9'j o p'dp'Qx^p'^^—z; (A.59) 
*FAx n (t,Z)=F- 1 [Q$ AXn (u>,iZ)], n = 0,l,2,..., (A.60) 

where A is the wave length and R — \x — x'\. represents the inverse Fourier 

transform. The first and second terms in eq.(A.59) represent respectively the con- 
tributions from high and low frequency components. We attached the symbol > to 
^fax„ meaning as before that the wave is superluminal. This is justified from the 
results of the experiment described in section 2. 

A6. The Donnelly-Ziolkowski Method (DZM) for Designing Subluminal, 
Luminal and Superluminal UPWs Solutions of the HWE and the Klein- 
Gordon Equation (KGE) 

Consider first the HWE for $ [eq.(A.2')j in a homogeneous medium. Let <&(a>, k) 
be the Fourier transform of <£>(£, x), i.e., 

~ -» r r+oo 

&(w,k)= d 3 x dt^{t,x)e- i{ ^-^ , (A.61a) 

JR 3 J-oo 
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$(t,x) = 7 - TI d 3 k du$(u,k)e i ( ki - ut \ (A.61b) 

(27TJ 4 JR 3 J-oo 

Inserting (A. 6 la) in the HWE we get 

(cu 2 -k 2 )$(u,k) = (A.62) 

and we are going to look for solutions of the HWE and eq.(A.62) in the sense of 
distributions. We rewrite eq.(A.62) as 

{u 2 - k\ - tt 2 )${u, k) = 0. (A.63) 

— * 

It is then obvious that any $(u;, k) of the form 

fc) = S(fi, /3) 5[cu - (/3 + fi 2 /4/3)] 5[A; Z -((3- n 2 /A(3)} , (A.64) 

where E(fl, (3) is an arbitrary weighting function, is a solution of eq.(A.63) since the 
5-functions imply that 

cu 2 - k 2 z = VI 2 . (A.65) 

In 1985' 3 °] Ziolkowski found a luminal solution of the HWE called the Focus 
Wave Mode. To obtain this solution we choose, e.g., 

E FWM (n,(3) = ^exp(-fi 2 V4/3), (A.66) 
whence we get, assuming (3 > and z > 0, 

Despite the velocities v\ = +1 and v 2 = —1 appearing in the phase, the modulation 
function of &fwm has very interesting properties, as discussed in details in ' 46 1. It 
remains to observe that eq.(A.67) is a special case of Brittingham's formula J 26 ' 

Returning to eq.(A.64) we see that the 5-functions make any function of the 
Fourier transform variables u, k z and Q to lie in a line on the surface u 2 — k 2 z — Q 2 = 
[eq.(A.63)j. Then, the support of the (^-functions is the line 

00 = (3 + Q 2 /4(3; k z = (3- Q 2 /4(3 . {AM) 

The projection of this line in the (cu, k z ) plane is a straight line of slope —1 ending 
at the point ((3,(3). When (3 = we must have fl — 0, and in this case the line 
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is uj = k z and x) is simply a superposition of plane waves, each one having 
frequency uj and traveling with speed c = 1 in the positive z direction. 

Luminal UPWs solutions can be easily constructed by the ZMt 46 ], but will not 
be discussed here. Instead, we now show how to use ZM to construct subluminal 
and superluminal solutions of the HWE. 

First Example: Reconstruction of the subluminal Bessel Beams $j o and the su- 
perluminal &xbb (AT- wave) 

Starting from the "dispersion relation" uj 2 — k 2 z — VL 2 = 0, we define 

$(cu,k) =Z(k~,r])5(k z -k cos r])5(uj-k). (A.69) 

This implies that 

k z = kcosrj; cosr/ = k z /uj, uj > 0, — l<cos?7<l. (A. 70) 

We take moreover 



Q = ksini]; k > . 



(A.71) 



— * — * 

We recall that Q = (k x ,k y ), p = (x,y) and we choose fl.p = flp cos 9. Now, 
putting eq.(A.69) in eq.(A.61b) we get 



\ roo 

(*, x) = -—- dk k sin 2 rj 



I-2TT 

/ d6E(k,r])< 
Jo 



,ikp sin Tj cos 



i(k cosr/z—kt) 



Choosing 



- p -kzQ sin»7 

Z(k,v) = (2-f-^ r - > 

k sin 77 



(A.72) 



(A.73) 



where z > is a constant, we obtain 

roo — 

$(t,x) = z sin r] / dke- kzosini > 
Jo 

Calling zosini] = a > 0, the last equation becomes 



1 r 2w 
2tt Jo 



a ik(cosr] z—t) 



poo — — 

$> o (t,f) = Ooj dke- kai \J (kpsmr ] )e lk(cosr ' z ~ t) . 
Writing k = k < and taking into account eq.(A.41) we see that 

J (£<psin?7)e^< (zcos, '-* ) 



(A.74) 



(A.75) 
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is a subluminal Bessel beam, a solution of the HWE moving in the positive z di- 
rection. Moreover, a comparison of eq.(A.74) with eq.(A.52) shows that (A. 74) is 
a particular superluminal X-wave, with B(k<) = e~ a ° k< . In fact it is the &xbb 
UPW given by eq.(A.56). 

Second Example: Choosing in (A. 72) 

Z(k,ri) = (27r) 3 e-^l COST '^ cot ry (A.76) 



gives 



POO _ _ 

$ > (t,x) = cos 2 77 / dkke- Z0 \ cos ^ k jQ{kpsmr])e- ik{ - cos ^ z -^ (A.77a) 
Jo 

[zq — isgn(cosi])(z — t/ cos rj)} 



[p 2 tan 2 77 + [zq + ?sgn(cosr7)(z — i/ cos if) ] 2 ] 3 / 2 ' 



(A. 77b) 



Comparing eq.(A.77a) with eq.(A.52) we discover that the ZM produced in this 
example a more general wave where B(k < ) = e - z o\ cos v\k< _ Obviously a?) 
given by eq.(A.77b) moves with superluminal speed (I/costi) in the positive or 
negative z-direction depending on the sign of cos 77, denoted sgn(cos?7). 

In both examples studied above we see that the projection of the supporting 
line of eq.(A.69) in the (u,k z ) plane is the straight line k z /u = cost], and cos 77 is 
its reciprocal slope. This line is inside the "light cone" in the (u, k z ) plane. 

Third Example: Consider two arbitrary lines with the same reciprocal slope that 
we denote by v > 1, both running between the lines u = ±k z in the upper half plane 
uo > and each one cutting the cu-axis at different values (3± and (3 2 (Fig. (12)). The 
two lines are projections of members of a family of HWE solution lines and each one 
can be represented as a portion of the straight lines (between the lines u = ±k z ) 

k z = v(lu - ft), k z = v(u - (3 2 ). (A.78) 

It is clear that on the solution line of the HWE, Q takes values from zero up to a 
maximum value that depends on v and (3 and then back to zero. 

We see also that the maximum value of f2, given by j3v / Vv 2 — 1, on any HWE 
solution line occurs for those values of u and k z where the corresponding projection 
lines cut the line u = vk z . It is clear that there are two points on any HWE solution 
line with the same value of Q in the interval 

< n < v(3/Vv 2 - 1 = n . (A.79) 
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It follows that in this case the HWE solution line breaks into two segments, as is the 
case of the projection lines. We can then associate two different weighting functions, 
one for each segment. We write 



{v 2 - 1) 



[f3v 2 + ^/v 2 (3 2 - Q 2 {v 2 - 1)] 



k z 



v\P+ ^f(3 2 v 2 -tt 2 {v 2 -1) 



X 



xS 



UJ — 



+ 



+ E 2 (n,v,(3)S 



(v 2 - 1) 

v[(3- ^(3 2 v 2 -tt 2 {v 2 - 1)] 
(v 2 - 1) 



X 



x 5 < uj - 



[pv 2 - ^Jv 2 p 2 -n 2 (v 2 -l)] 

(v 2 - 1) 



(A.80) 



Now, choosing 



Si(ft,u,/3) = Z 2 (n,v,P) = (2tt) 3 /2^ _ Q 2 



we get 
Then 

$ v ,n(t,p, z) = exp 



'ij3v{z — vt) 
v(z - vt) 



dX XJ Q (Q pX) cos 



Qov (z — t/v) 



id 



{^^{z-t/vf + p 2 } 



(A.81) 



If we call v< — - < 1 and taking into account the value of VLq given by eq.(A.79), 

v 

we can write eq.(A.81) as 



(t,P,z) 



sin(fi ^<) a 



iflo(z-vt) . 



nl/2 



a: 2 + y 2 + 



l-vl 



(z- V< tf 



(A.82) 



which we recognize as the subluminal spherical Bessel beam of section A4 [eq.(A.31)j. 
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Figure 12: Projection of the support lines of the transforms of two members of a 
family of subluminal solutions of the HWE. 

Klein-Gordon Equation (KGE): We show here the existence of subluminal, lu- 
minal and superluminal UPW solutions of the KGE. We want to solve 

^ - V 2 + m 2 ^j <S> KG {t,x) = 0, m>0, (A.83) 

with the Fourier transform method. We obtain for & KG (u, k) (a generalized func- 
tion) the equation 

{uj 2 - k 2 z - (Q 2 + m 2 )}$ KG (u, k) = 0. (A.84) 

As in the case of the HWE, any solution of the KGE will have a transform 
$(u,k) such that its support line lies on the surface 

lo 2 - k\ - (n 2 + m 2 ) = . (A.85) 

From eq.(A.85), calling f2 2 + m 2 = K 2 , we see that we are in a situation identical 
to the HWE for which we showed the existence of subluminal, superluminal and 
luminal solutions. We write down as examples one solution of each kind. 
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Subluminal UPW solution of the KGE. To obtain this solution it is enough to 



change in eq.(A.81) VL = vfi/^v^l -> tt£ G 



vf3 



^ G (t,p,z) = exp 



i(3v(z - vt) \ sin(£ffi G £<) 



m 



1/2 



. We have 



1 



+ r + - Az-v<t) 



l-v% 



1/2 



, v< = l/v. (A.86) 



Luminal UPW solution of the KGE. To obtain a solution of this type it is enough, 
as in eq.(A.64), to write 



<b KG = E(Q, f3)8[k z - (VI 2 + (m 2 - f3 2 )/2f3)]8[u - (fi 2 + (m 2 + (3 2 )/2(3)] . (A.87) 
Choosing 

E(Sl,/3) = ^-exp(-z n 2 /2P), z > 0, (A.88) 

gives 

*?°(*,*) = 

PY r,/_n 2 /^/9rr„ _ ,7r - ^11 

-. (A.89) 



= exp(^(m 2 - (3 2 )/2(3) exp(-tt(m 2 + (3)/2(3)^ { p ^/ 2 ^ ^ ^ 



[zq -i(z- t)} 

Superluminal UPW solution of the KGE. To obtain a solution of this kind we 
introduce a parameter v such that < v < 1 and write for <& KG in (A. 84) 



$« G (tu,n,k z ) = E(n,v,p)s 



-pv 2 + ^J(tt 2 + m 2 )(l-v 2 )+v 2 (3 2 ) 



x5 



1-v 2 

v + ^J(tt 2 + m 2 )(l-v 2 )+v 2 [3 2 ) 



x 



(A.90) 



Next we choose 



(2n) 3 exp(-z ^n 2 + Q 2 ) 

^n 2 + n 2 



(A.91) 
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where zq > is an arbitrary parameter, and where 



2„,2 



n 2 = pL- + m 2 . (A.92) 



1 - u 



Then introducing t>> = l/t> > 1 and 7> = , we get 

v ?; > - 1 



S) k G>( . ^ _ pvn f - ^ 2 )(^ - 1 g^goyg - *r>(* - y>t)} 2 + + I/ 2 } 

?j Q \ L : - 1 ) — ^Xp \ T, ( / j 

,P I P v J /[ Z „-i7>(2-«>f + ^ + !/ 2 

(A.93) 

which is a superluminal UPW solution of the KGE moving with speed v > in the z 
direction. From eq.(A.93) it is an easy task to reproduce the superluminal spherical 
Bessel beam which is solution of the HWE [eq.(A.30)j. 

A7. On the Energy of the UPWs Solutions of the HWE 

Let <S> r (t, x) be a real solution of the HWE. Then, as it is well known, the energy 
of the solution is given by 



dv 

M 3 



+ lim / / dS$ r n.V$ r , (A.94) 

R^ooJ Js(R) 



where S(R) is the 2-sphere of radius R . 

We can easily verify that the real or imaginary parts of all UPWs solutions of the 
HWE presented above have infinite energy. The question arises of how to project 
superluminal waves, solutions of the HWE, with finite energy. This can be done if 
we recall that all UPWs discussed above can be indexed by at least one parameter 
that here we call a. Then, calling <& a (t,x) the real or imaginary parts of a given 
UPW solution we may form "packets" of these solutions as 

= j daF(a)$ a (t,x) (A.95) 

We now may test for a given solution $ Q and for a weighting function F(a) if 
the integral in eq.(A.94) is convergent. We can explicitly show for some (but not 
all) of the solutions showed above (subluminal, luminal and superluminal) that for 
weighting functions satisfying certain integrability conditions the energy e results 
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finite. It is particularly important in this context to quote that the finite aperture 
approximations for all UPWs have, of course, finite energy. For the case in which $ 
given by eq.(A.95) is used to generate solutions for, e.g., Maxwell or Dirac fields (see 
Appendix B), the conditions for the energy of these fields to be finite will in general 
be different from the condition that gives for <3> a finite energy. This problem will 
be discussed with more details in another paper. 

Appendix B. A Unified Theory for Construction of UPWs 
Solutions of Maxwell, Dirac and Weyl Equations 

In this appendix we briefly recall the main results concerning the theory of Clif- 
ford algebras (and bundles) and their relationship with the Grassmann algebras 
(and bundles). Also the concept of Dirac- Hestenes spinors and their relationship 
with the usual Dirac spinors used by physicists is clarified. We introduce moreover 
the concepts of the Clifford and Spin-Clifford bundles of spacetime and the Clifford 
calculus. As we shall see, this formalism provides a unified theory for the construc- 
tion of UPWs subluminal, luminal and superluminal solutions of Maxwell, Dirac and 
Weyl equations. More details on the topics of this appendix can be found in' 6 " 9 ! 
and[ 81 l 

Bl. Mathematical Preliminaries 

Let M. = (M,g,D) be Minkowski spacetime. (M,g) is a four-dimensional 
time oriented and space oriented Lorentzian manifold, with M ~ R 4 and g e 
sec(T*M x T*M) being a Lorentzian metric of signature (1,3). T*M [TM] is 
the cotangent [tangent] bundle. T*M = U xeM T*M and TM = U xeM T x M, and 
T X M ~ T*M ~ J? 1 ' 3 , where J? 1 ' 3 is the Minkowski vector space' 60 ' 61 - 62 !. £> [ s the 
Levi-Civita connection of g, i.e., Dg = 0, T(D) = 0. Also R(D) = 0, T and R 
being respectively the torsion and curvature tensors. Now, the Clifford bundle of 
differential forms C£(M) is the bundle of algebras Cl{M) = U xeM Ci(T*M), where 

G M, Ci(T*M) = a h3 , the so called spacetime algebra' 9 ' 81 " 851 . Locally linear 
space over the real field R, Ci(T*(M)) is isomorphic to the Cartan algebra /\(T*M) 
of the cotangent space and A{T*M) = £Lo A k {T*M), where f\ k {T*M) is the Q- 
dimensional space of /c-forms. The Cartan bundle /\(M) = U xE m A(T*M) can then 
be thought as "embedded" in Cl{M). In this way sections of Ct(M) can be repre- 
sented as a sum of inhomogeneous differential forms. Let {e M = ^-} e sec TM, 
(/i = 0,1,2,3) be an orthonormal basis g(e^,e u ) = rj^ = diag(l, — 1, — 1, — 1) and 
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let {y = dx v } G sec h}(M) C sec C£(M) be the dual basis. Then, the fundamental 
Clifford product (denoted in what follows by juxtaposition of symbols) is generated 
by YY + YY = 2^ and if C G secG?(M) we have 

c = s + + ^yy + ^ p yyy + > (b.i) 

where Y = 7°7 1 7 2 7 3 = dx°dx l dx 2 dx 3 is the volume element and s, u„, b 



fll U [IV) "'ftupt 



p G sec A (M) C sec (X(M). For A r G secA r (M) C sec C£(M), B s G secA s (M) 
we define' 9,82 ' A r ■ B s = (A r B s )\ r ^ s \ and A r A B s = (A r B s ) r+s , where ( )k is the 
component in /\ h (M) of the Clifford field. 

Besides the vector bundle C£(M) we also need to introduce another vector bundle 
<^s P m + (i,3)(M) [Spin + (1,3) ~ SL(2,(T)] called the Spin-Clifford bundle! 8 - 81 ' 84 ]. We 
can show that G?s P in + (i,3)(M) ~ C£(M)/1Z, i.e. it is a quotient bundle. This means 
that sections of C£spin + (i,3)(M) are equivalence classes of sections of the Clifford 
bundle, i.e., they are equivalence sections of non-homogeneous differential forms 
(see eqs. (|B~2l , |B~3D below). 

Now, as is well known, an electromagnetic field is represented by F G sec /\ 2 (M) C 
sec Ci(M). How to represent the Dirac spinor fields in this formalism ? We can show 
that the even sections of C£s P m + (i,3)(M), called Dirac-Hestenes spinor fields, do the 

job. If we fix two orthonormal basis £ = {7^} as before, and £ = {7^ = Rj^R = 
A£Y} with A£ G SO+(l,3) and R(x) G Spin + (1,3) Vx G M, RR = RR = 1, and 
where ~ is the reversion operator in Gt\^, then! 8 ' 81 ! the representatives of an even 
section if) G secCfspin + (i,3)(^0 are the sections ip^ and ip^ of Cl{M) related by 

^ = ^R (B.2) 

and 

^ = s + ^b llu YY + Pl 5 - (B.3) 

Note that has the correct number of degrees of freedom in order to represent a 
Dirac spinor field, which is not the case with the so called Dirac-Kahler spinor field 

(seel 8 ' 81 !). 

Let * be the Hodge star operator * : A k (M) -> A 4 ~ k (M). We can show that 
if A p G secA p (-^0 C secG?(M) we have *A = AY- Let d and 5 be respectively 
the differential and Hodge codifferential operators acting on sections of f\(M). If 
ujp G sec A P (M) C secd(M), then 5lu p = (-l) p d*co> p , with = identity. 

The Dirac operator acting on sections of Ci(M) is the invariant first order dif- 
ferential operator 

9 = YD,, (B.4) 
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and we can show the very important result (see e.g.' 6 '): 

d = dA+d- = d-5. (B.5) 
With these preliminaries we can write Maxwell and Dirac equations as follows' 82 ' 85 ': 

OF = 0, (B.6) 

d^ s7 V + n#E7° = 0. (B.7) 

We discuss more this last equation (Dirac-Hestenes equation) in section B.4. If 
m = we have the massless Dirac equation 

d^E = 0, (B.8) 

which is Weyl's equation when ip^ is reduced to a Weyl spinor field (see eq. ([B.12| be- 
low). Note that in this formalism Maxwell equations condensed in a single equation! 
Also, the specification of ips depends on the frame S. When no confusion arises we 
represent simply by ip. 

When V'sV's 7^ 0, where ~ is the reversion operator, we can show that has 
the following canonical decomposition: 

^ = ^pe^ 2 R, (B.9) 

where p, (3 G secA°(M) C secC^(M) and R G Spin + (1,3) C Ci\^ Vi 6 M. (3 is 
called the Takabayasi angle ®. 

If we want to work in terms of the usual spinor field formalism, we can translate 
our results by choosing, for example, the standard matrix representation of {7 M }, and 
for given by eq. ( |B.3[ ) we have the following (standard) matrix representation' 8 ' 49 ] : 



(B10) 



where 



s - ib 12 &13 - ib 2 3 \ , _ ( -b 3 + W -&oi + *&oa 
-613 - ib 23 s + ib 12 ) ' 2 I -b i - ib m b 03 + ip 



(B.11) 



with s, &12, • • • real functions; * denotes the complex conjugation. Right multiplica- 
tion by 

( 1 \ 




V o / 
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gives the usual Dirac spinor field. 

We need also the concept of Weyl spinors. By definition, ip G sec C£ + (M) is a 
Weyl spinor iff 83 ' 

7 5 ^ = ±^721 • (B.12) 

The positive [negative] "eigenstate" of 75 will be denoted ip + For a general 

ip G sec Ci + (M) we can verify that 

^± = \[^T 75^721] (B.13) 

are Weyl spinors with eigenvalues ±1 of eq. flB.12p . 

We recall that the even subbundle Cl + (M) of C£(M) is such that its typical fiber 
is the Pauli algebra C£ 3>0 = Cl\^ (which is isomorphic to <F(2), the algebra of 2 x 2 
complex matrices). The isomorphism Ci 3 ,o = Cif 3 is exhibited by putting <jj = 7^70, 
whence o~i<jj + ofti = 25^. We recall also' 8 ' 81 ' that the Dirac algebra is Cl^\ = <T(4) 
(see section B4) and fl? 4 ,i = <F ® <^i,3 [86] - 

B2. Inertial Reference Frames (/), Observers and Naturally Adapted 
Coordinate Systems 

Let Al = (M,g,D) be Minkowski spacetime. An inertial reference frame (irf) 
J is a timelike vector field I G sec TM pointing into the future such that g(I, I) = 1 
and DI = 0. Each integral line of I is called an inertial observer. The coordinate 
functions (x^), /x = 0, 1, 2, 3 of the maximal atlas of M are said to be a naturally 
adapted coordinate system to / (nacs/J) if / = d/dx° [ 61>62 1. Putting / = cq we can 
find ej = d/dx\i = 1,2,3 such that g(e^,e u ) = r]^ u and the coordinate functions 
x^ 1 are the usual Einstein-Lorentz ones and have a precise operational meaning: 
x° = ct^*\ where t is measured by "ideal clocks" at rest on I and synchronized "a la 
Einstein", x\i = 1,2,3 are determined with ideal rules' 61 ' 62 ]. (We use units where 
c = l.) 

B3. Maxwell Theory in Cl{M) and the Hertz Potential 

Let e M G sec TM be an orthonormal basis, g(e fl ,e u ) = rj^ and e M = d/dx 11 
(fi,u = Q,l, 2, 3), such that e determines an IRF. Let 7^ G sec A 2 (M) C sec Cl(M) 

(*'c is the constant called velocity of light in vacuum. In view of the superluminal and sub- 
luminal solutions of Maxwell equations found in this paper we don't think the terminology to be 
still satisfactory. 
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be the dual basis and let 7 M = i] IIU '-f' / be the reciprocal basis to 7 M , i.e., 7^.7^ = 8%. 
We have 7 M = dx^ . 

As is well known the electromagnetic field is represented by a two-form F G 
secA 2 (M) C sec C£(M). We have 



F = -F^7„7„ F>» 



( 


—E 1 


—E 2 


—E 3 \ 


E 1 





-B 3 


B 2 


E 2 


B 3 





-B 1 


{ E 3 


-B 2 


B 1 


J 



;b.i4) 



where (E l ,E 2 ,E 3 ) and 5 2 , i? 3 ) are respectively the Cartesian components of 
the electric and magnetic fields. Let J G secA 1 (^) C sec C£(M) be such that 



J = ^7/, = P7o + J l li + A2 + J 3 7s, 



(B.15) 



where p and (J 1 , J 2 , J 3 ) are respectively the Cartesian components of the charge 
and of the three-dimensional current densities. 

We now write Maxwell equation given by |B.6| in C£ + (M), the even sub-algebra 
of C£(M). The typical fiber of C£ + (M), which is a vector bundle, is isomorphic to 
the Pauli algebra (see section Bl). We put 



0* = 7*7o, 1 = o"iO"20"3 = 7o7i7273 = 7s- 



(B.16) 



Recall that i commutes with bivectors and since i 2 = — 1 it acts like the imagi- 
nary unit % = y— T in C£ + (M). From eq.( |B.14|) , we get 



F = E + \B 



;B.17) 



with E = E l di, B = B j tjj, i,j = 1,2,3. Now, since d = 7^(9^ we get d"fo 
d/dx° + did 1 = d/dx° — V. Multiplying eq. (|B.6|) on the right by 70 we have 



97070-F70 = Jlo, 

(d/dx° - V)(-fi + iB) = p + J, 
where we used 7°fi7o = —E + iB and J = J l Bi. From eq.( p.l8|) we have 



(B.18) 



-d E + id B + VE - iVfi = p + J (B.19) 
-d E + id B + V.fi + V A E - iV.fi -iVAfi = p + J (B.20) 
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We have also 

-iVAA = Vxl (B.21) 

since the usual vector product between two vectors a = a % di, b = h l a,i can be 
identified with the dual of the bivector a A b through the formula a x b = — i(a A 
b). Observe that in this formalism a x b is a true vector and not the meaningless 
pseudovector of the Gibbs vector calculus. Using eq. QB.21J) and equating the terms 
with the same grade we obtain 

V.E = p; V x B -d E = J; 

(B.22) 

Vx£ + d B = ; V .B = ; 

which are Maxwell equations in the usual vector notation. 

We now introduce the concept of Hertz potential^ 19 ' which permits us to find 
nontrivial solutions of the free "vacuum" Maxwell equation 

OF = (B.23) 

once we know nontrivial solutions of the scalar wave equation, 

□$ = (d 2 /dt 2 - V 2 )$ = 0; $ G sec A°(M) C sec C£(M) . (B.24) 

Let A G sec/^^M) C sec C£(M) be the vector potential. We fix the Lorentz 
gauge, i.e., d.A = —5A = such that F = OA — (d — 5) A = dA. We have the 
following important result: 

Theorem: Let n G secA 2 (^0 C secC£(M) be the so called Hertz potential. If ir 
satisfies the wave equation, i.e., Ott = d 2 7t = (d — 5)(d — 5)n = —(d5 + 5d) tt = 
and if we take A = —5tt, then F = dA satisfies the Maxwell equation OF = 0. 

The proof is trivial. Indeed, A = —5tt, implies 5 A = —5 2 n = and F = dA = 
dA. Then OF = (d — 5)(d — 5) A = 5d{5'K) = —5 2 dn = 0, since 5dir = —d5ir from 
d 2 ir = 0. 

From this result we see that if $ G secA°(^0 C sec C£(M) satisfies <9 2( £ = 0, 
then we can find non trivial solution of dF = 0, using a Hertz potential given, e.g., 
by 

vr = $ 7l72 . (B.25) 

In section 3 this equation is used to generate the superluminal electromagnetic X- 
wave. 
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We now express the Hertz potential and its relation with the E and B fields, in 
order for our reader to see more familiar formulas. We write 7r as sum of electric 
and magnetic parts, i.e., 



7T = 7T e + \lt. n 



7T e = -7T 0i <7;, 7T m = -TT 23 ^ + TT 13 ^ - 7T 12 <7 3 



Then, since A = dir we have 



E e = V x V x 7? e ; 5 e = -9 (Vx 7r e ) ; 
E m = -<9 (V x 7? m ) ; 4=-VxVx vr r , 



'B.26) 



A = ^(dn-nd) (B.27) 
A 7o = -<9 7f e + V.vf e - (V x 7? m ) (B.28) 

and since A = A M 7 M we also have 

A° = V.7f e ; A=A% = -^7f e - V X 7? m . 

Since £ = - VA° - 7^77^ and B = VxI,we obtain 
ox" 

E = -5b(V x 7f m ) + V x V x 7f e ; (B.29) 

B = VX (-9 7T e - V X 7T m ) = -5b(V X 7T e ) - V X V X 7T m . (B.30) 

We define E e ,B e ,E m , B m by 



(B.31) 



We now introduce the 1-forms of stress-energy. Since dF = we have Fd = 0. 
Multiplying the first of these equation on the left by F and the second on the right 
by F and summing we have: 

{\/2){FdF + FdF) = d^{l/2)F^F) = = 0, (B.32) 

where Fd = -(d^\F a ^ a ^)^ . Now, 

- \{F^F)Y = -\{F^FY) (B.33) 
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Since ^.F = -(^F - F^) = F.^, we have 

V" = -((F.r)FY)o-l(YF 2 Y)o 

= -(^.(F.yj-^F.Fjy.y (b.34) 

= F^F^ VaX + ^rTF aP F^, 

which we recognize as the stress-energy momentum tensor of the electromagnetic 
field, and T M = T^ v lv _. 

By writing F = E + iB as before we can immediately verify that 

T = -^F l0 F 



-(E 2 + B 2 ) + (E x B) 
2 



7o- (B.35) 

We have already shown that <9 M T M = 0, and we can easily show that 

d.T" = 0. (B.36) 

We now define the density of angular momentum. Choose as before a Lorentzian 
chart (x 11 ) of the maximal atlas of M and consider the 1-form x = x^^^ = x^^. 
Define 

M„ = x A = ^(x a T^ - x u T a ^ a A 7 " ■ 
It is trivial to verify that as T^ v = T Vil and d^T^ = 0, it holds 

d»M„ = 0. (B.37) 
The invariants of the electromagnetic field F are F. F and F A F and 

F 2 = F. F + F A F; 

F. F — —F^F^ ; F A F = -^F^F^e^p . (B.38) 
Writing as before F = E + iB we have 

F 2 = (£ 2 - 5 2 ) + 2iE.il = F. F + F A F. (B.39) 
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B4. Dirac Theory in C£(M) 

Let S = {7 M } G secA 1 !^) C sec C£(M) be an orthonormal basis. Let ^ G 
sec(A°(M) + A 2 (M) + A 4 (M)) C sec C£(M) be the representative of a Dirac- Hestenes 
Spinor field in the basis E. Then, the representative of Dirac equation in C£(M) is 
the following equation (H = c = 1): 

<9V>s7i72 + mipxjo = . (B.40) 

The proof is as follows: 

Consider the complexification C£c{M) of C£(M) called the complex Clifford bun- 
dle. Then C£c{M) = (F <S> C£(M) and by the results of section Bl it is trivial to see 
that the typical fiber of C£c(M) is Cl^\ = <T <g> C^i,3, the Dirac algebra. Now let 
{r , r 1; T 2 , T 3 , r 4 } C secA 1 ( j ^) C sec CEc(M) be an orthonormal basis with 

r a r b + r 6 r a = 2^ , (B.4i) 

5f afe = diag(+l, +1, +1, +1, -1) . 

Let us identify 7 M = r At r 4 and call / = ror 1 r 2 r 3 r4. Since I 2 = — 1 and 7 
commutes with all elements of C£^i we identify I with i = 1 and 7^ with the 
fundamental set of C£(M). Then if ^4 G secC£c(^0 we have 

.A = $ s + A£ T/i + ^^7m% + ^'7^7, + $ pT5, (B.42) 

where $ s , $ p , A£, ££7 r^ p G sec d? <g> A°(M) C sec C£ C (M), i.e., Vx G M, $ s (z), 
B(?(x), r^ up (x) are complex numbers. 

Now, 

/=i(i + 70)^(1 + ^7172); / 2 = /, 

is a primitive idempotent field of C£q{M). We can show that if = 7271/. From 
( [B.40Q we can write the following equation in C£c(M): 

<9^s7i72/ + mipslof = (B.43) 
difeif - mifcf = (B.44) 

and we have the following equation for \l/ = ipY,f: 

id^ - m* = 0. (B.45) 
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Using for 7„ the standard matrix representation (denoted here by 7 ) we get 

— M 

that the matrix representation of eq. ( B.45|) is 



i^fd^) - m|*> = (B.46) 

where now |$) is a usual Dirac spinor field. 

We now define a potential for the Dirac- Hestenes field ipv Since ?/>s G sec C£ + (M) 
it is clear that there exist A and -B G secA C sec C£(M) such that 

Ve = <9(A + 7 5j B), (B.47) 



since 



d(A + lb B) = d.A + dAA — ^d.B - A B (B.48) 

= S + B + l5 P- (B.49) 

5 = B = OA A- l5 dAB; P = -d.B. 

We see that when m = 0, ipY, satisfies the VFey/ equation 

dipx = . (B.50) 



Using eq.( |B.50| ) we see that 

<9 2 A = <9 2 £ = 0. (B.51) 

This last equation allows us to find UPWs solutions for the Weyl equation once 
we know UPWs solutions of the scalar wave equation □$ = 0, $ G secA°(^0 C 
sec C£(M). Indeed it is enough to put A = (A + 75-B) = $(1 + 75)17, where v is a 
constant 1-form field. This result has been be used in' 48 ] to present subluminal and 
superluminal solutions of the Weyl equation. 

We know (see appendix A, section A5) that the Klein-Gordon equation have 
superluminal solutions. Let $> be a superluminal solution of □$> + m 2 $> = 
0. Suppose $> is a section of C£c(M). Then in C£c{M) we have the following 
factorization: 

(□ + m 2 )$ = (d + im)(d - im)$ = 0. (B.52) 

Now 

V > = (d-im)$ > f (B.53) 
is a Dirac spinor field in C£c{M), since 

(d + im)^ > = (B.54) 
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If we use for $ in eg. ( |B.52|) a subluminal or a luminal UPW solution and then use 
eq. (|B.53|) we see that Dirac equation also has UPWs solutions with arbitrary speed 
< v < oo. 
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